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ABSTRACT 

The recent trend of using fine water mist systems to replace the legacy HALON- 
1301 fire suppression systems warrants further study into other applications of the water 
mist systems. Preliminary research and investigation indicates that fine mists (20-25 pm 
droplet size) may reduce peak overpressures of a shock wave traveling through a space. 
Such pressure reductions could be used to mitigate the destructive effects of a shock wave 
(initiated by an explosive device) traveling through a structure. 

Currently these blast mitigation effects have only been demonstrated in small-scale 
shock tube tests and computer simulations. Uncertainty exists as to the scalability of such a 
system. The intention of this research is to investigate the applicability of such a blast 
mitigation system for shipboard use. Study into the degree of mitigation necessary to make 
a system practical for shipboard installation was conducted. In addition, a theoretical study 
of the mechanisms of blast mitigation using water mists was completed. 

Preliminary design of a full-scale system was examined. Given the recent trend 
toward tumblehome hull forms in future Naval Combatant designs, there exists strong 
applicability of this system in the “dead 5 " spaces created by the shaping of the tumblehome 
hull. Further work is needed in numerical modeling and laboratory testing of specific 
phases of the mitigation. The end goal is a feasible design of a blast mitigation system to be 
used in the outermost spaces of Naval Combatants to protect interior vital system spaces. 
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NOMENCLATURE 

DDS: U.S. Navy Design Data Sheets. Contain construction and performance 
requirements and guidelines for U.S. Navy ships. 

GenSpecs: U.S. Navy General Specifications for construction and performance of 
U.S. Navy ships. 

Survivability: Concept that includes all aspects of protecting personnel, weapons, 
and supplies while simultaneously deceiving the enemy. 

ODS: Ozone Depleting Substance. 

ODP: Ozone Depletion Potential. 

CFC: Chlorofluorocarbon, a common refrigerant. 

CPSMA: Commission on Physical Sciences, Mathematics, and Applications. 

T v : Particle equilibration time, 
m: Individual particle mass. 

T|: mass loading ratio. 

P D : Peak incident overpressure. Peak pressure in the incident shock wave. 

P a : Ambient atmospheric pressure. 

P r : Maximum reflected overpressure. Peak pressure due to additive effect of 
incident and reflected shock waves. 

q Q : Dynamic pressure. Accounts for dynamic effects of pressure pulse duration. 

P s : Stagnation pressure. Static equivalent pressure including dynamic effects. 
x: Duration of positive phase of pressure pulse. 

Z: Scaled blast range. 

W: Weight of explosive for blast calculations. 








|Li: Ductility factor. Used in DDS 100-9 to account for allowable plastic 
deformation levels. 

Rpi Panel resistance in pounds force. Used in DDS 100-9 to evaluate stiffened 
panel design. 

a G : Plastic flow stress 

M 0 : Plastic collapse moment 

CT odyn : dynamic yield stress 

8 dot : strain rate (= rupture strain/pulse time) 

L: Frame spacing. Horizontal length of stiffened panel. 

H: Plate thickness. 

B: Stiffener spacing or deck height. 

PVLS: Peripheral Vertical Launch System 




INTRODUCTION 


Central to the strength of a naval force is its ability to survive and continue to fight 
through a batde situation. The United States Department of Defense Dictionary of 
Military Terms defines survivability as “a concept that includes all aspects of protecting 
personnel, weapons, and supplies while simultaneously deceiving the enemy.” The ability 
of a system to retain functionality in a damaged state is central to the concept of 
survivability and therefore is a primary concern in the design of naval combatants. The 
protection of personnel and equipment from battle damage is essential to the effectiveness 
of a ship. Warships must be designed with inherent resistance to hull and shipboard system 
damage inflicted under batde conditions. Part of this protection includes methods of 
strengthening the structure of the ship to resist damage from blast-type loading. The 
stronger die ship, the more likely it is to survive batde damage and continue its ability to 
fight in a damaged condition. 

Strengthening ships to blast loading is traditionally achieved through enhanced 
structural design, but can include measures to mitigate the peak overpressures and loads 
placed on a ship by shock impingement. Water mists have been shown to achieve a degree 
of mitigation of blast shocks in a confined space, namely a reduction in the velocity and 
peak overpressure of the impinging shock wave. The potential application of this feature of 
water mists to U.S. Naval combatants is presented in this study. 

Ship Design for Blast Resistance 

There are obvious trade-offs to batde hardening of Naval assets. Allocating 
significant amounts of weight and space to armor plating or blast-resistant structures will 
cost in speed, efficiency, endurance, and overall effectiveness of the platform. In the days 
before technologically advanced radars and guided weapons, armor plating was essential to 
ship design since major threats included ship-to-ship gunfire and close-in ballistic fire. 
However, recent years have turned the focus on creating ships with low signatures and 
small radar cross sections to make targeting by enemies more difficult. To accomplish 




these goals, ships must appear small and quiet, often conflicting design features to side-hull 
protection. Modem warships are also filled with the sensitive electronic equipment 
required to compete in today’s high-tech arenas. These ships contain a much higher ratio 
of combat system equipment weight to ship structural weight in comparison to naval 
combatants of the past. The next generation of ships is entirely electric and must carry a 
previously unheard of amount of electrical generation and distribution equipment to 
complete their required missions. They will operate in an environment where the force that 
can conduct the best surveillance, advanced warning, and precision strike missions will 
dominate the field. This means that ships must be fast, technologically advanced, and 
robust. 

As technologically advanced next-generation ships are designed to meet the 
computing and combat systems needs of the future, the concept of survivability in ship 
design has changed. Battle-hardening no longer means the ability to withstand hits from 
explosive shells. Ships are designed to project power far away and detect and combat 
threats from equally long distances. As a result, heavy armor has been traded for high-tech 
detection and engagement systems. However, the close-in threat has not disappeared. The 
most daunting threats to the new generation of ships are from undetected, fast moving 
cruise missiles and unexpected terrorist plots. Ships are designed to keep a full operational 
picture of an entire battle space, so the threat of direct fire from another surface combatant 
or strike from an air asset is minimized. The new dangers to these ships are non-traditional 
war fighting scenarios such as small boat attacks when operating in port or in vulnerable 
situations close to shore. This close-in threat is evidenced by the attack on the USS COLE 
(DDG-67) in the port of Aden, Yemen on October 12,2000. 

USS COLE was performing a routine refueling operation in a protected harbor 

when a small vessel appearing to be a tender pulled alongside. An explosion detonated on 

the small vessel caused extensive damage to the ship. Unofficial estimates state that the 

explosive charge was equivalent to approximately 400-700 pounds of TNT. It was 

detonated along the waterline at a standoff distance of 0 to 10 feet from the hull. The 

explosion caused a 40 foot by 45 foot hole in the port side of the ship extending below the 
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waterline. Two longitudinal watertight subdivisions were breached, exposing the interior of 
the ship to fire and flooding damage well past centerline. Figure 1 is a photograph of the 
damage to the port side of the USS COLE, revealing the blast damage to the side hull of 
the ship. The resulting damage crippled the state-of-the-art Arleigh-Burke class destroyer 
and instigated a close look at methods of mitigating the damage caused by such attacks to 
Naval vessels. 



Figure 1: Port Side Damage to USS COLE (Photo from U.S. 
Navy Information Office) 


In a case such as the explosion alongside the USS COLE, the tearing of the hull 
and subsequent travel of the explosive shock wave through the ship infrastructure resulted 
in a complete loss of power and propulsion capabilities, along with significant injury and 
loss of life to personnel. A method of mitigating the spread of the shock wave through the 
ship might have isolated the damage to the outermost spaces adjacent to the blast, 
protecting the inboard engineering spaces. Preserving the integrity of these vital spaces 
would allow continued operation of propulsion and electrical generation engines and 
maintain war-fighting capabilities. 

Arleigh-Burke class destroyers such as USS COLE are not equipped with specific 
side-hull protection systems. They are designed to withstand the necessary ship loads in 
accordance with U.S. Navy Design Data Sheets (DDS). DDS’s are design guidelines that 
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outline a systems-based approach to the mitigation of damage for ship structures. They 
were developed using classified explosive deformation and failure studies for naval vessels, 
empirically based on data accrued through years of live fire tests conducted by the Office of 
Naval Research (ONR). The blast loading guidelines in these documents are primarily 
concerned with protecting the ship from the stresses produced by own-ship weapons 
systems and hull dynamic loading. For example, there are requirements to design 
bulkheads and decks to withstand the loading of the overpressure from a gun barrel firing 
or from the impingement of missile rockets as they launch. However, there are no specific 
requirements for U.S. Navy ships to have armor plating or blast protection to withstand 
damage from outside sources other than nuclear overpressures. Blast protection or 
mitigation capabilities, if installed, are merely design features of the specific platform. 

The next generations of U.S. Navy destroyers, DD(X), are designed with a bit 
more side-hull protection as a consequence of the location of their missile launchers. 
DD(X) has a tumblehome hull form; the sides of the hull slope inward from a maximum 
beam below the waterline. Vertical launch missile cells are positioned direcdy inboard of 
the side hull in this design, creating a dead space outboard of the missiles due to the slope 
of the hull. In order to protect the ship from damage due to unintentional detonation of 
the missiles housed in these launchers, current design uses reactive armor and heavy plating 
inboard of the launcher. This design is necessary to protect the ship, but cosdy in terms of 
space and weight. 

- An alternative means to partially reduce blast loading on the side hull of DD(X) 
could produce significant savings for the ship design. Blast mitigation would need to be 
performed with minimal impact on major ship systems, preferably using extensions of 
currendy installed systems. Water mists are currendy being used as fire suppression systems 
on the next-generation destroyer and future U.S. Naval combatants. It appears that some 
degree of mitigation may be achieved by proper application of water mist to a space prior 
to the propagation of a shock wave. Therefore, the feasibility of using water mist to 
mitigate blast damage must be investigated. 



HALON Replacement and Water Mist 

Current U.S. Naval combatants primarily rely on HALON-1301 as a space flooding 
fire protection system for vital shipboard spaces such as engine rooms. Halons are low- 
toxicity, chemically stable compounds that have been used for fire and explosion protection 
throughout this century. HALON-1301 (a gaseous agent) is used mainly in total flooding 
fire suppression systems. Halons have proven to be extremely effective fire suppressants, 
which are clean (leave no residue) and remarkably safe for human exposure. They work by 
breaking the chain reaction required for a fire to bum unabated. 

Three things must come together at the same time to start a fire. The first 
ingredient is fuel (anything that can bum), the second is oxygen (normal breathing air is 
ample) and the last is an ignition source (high heat can cause a fire even without a spark or 
open flame). Traditionally, to stop a fire you need to remove one side of the triangle - the 
ignition, the fuel or the oxygen. Halon adds a fourth dimension to fire fighting - breaking 
the chain reaction. It stops the fuel, the ignition and the oxygen from dancing together by 
chemically reacting with them. Many people believe that Halon displaces the air out of the 
area it is dispensed in, but this is false. Even for the toughest hazards, less than an 8% 
concentration of Halon by volume is required. There is still plenty of air in the space for 
personnel to use during the evacuation process. 

Under the Clean Air Act (CAA), the U.S. banned the production and import of 

virgin HALON-1301 beginning January 1, 1994 in compliance with the 1987 Montreal 

Protocol. The Montreal Protocol placed a ban on the production of ozone depleting 

substances (ODSs) for industrialized nations in January 1996. A compound’s ability to 

destroy ozone depends on many factors, including the amount of chlorine and/or bromine 

that it contains. To aid them in comparing compounds, scientists have developed a relative 

scale called the ozone depletion potential (ODP). Chlorofluorocarbons (CFCs), the 

common refrigerants found in refrigerators and car air conditioners, have been assigned the 

value 1 as a reference. HALON-1301 has the value between 10 and 16, meaning it has 10- 

16 times the more potential for destroying the ozone layer. Halon use worldwide is 

significantly less than that of CFCs, so even though it is more damaging to the ozone layer, 
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there is not as much of it released into the atmosphere. In fact, it is estimated that overall 
Halons account for less than 20% of ozone depletion. Recycled Halon and inventories 
produced before January 1,1994 are now the only sources of supply in the United States. 

There are a number of traditional fire extinguishing agents, such as water, carbon 
dioxide, dry chemicals, and foam that are good alternatives to Halons for many 
applications. In addition, recent research has led to the commercialization of new agents 
and technologies. These fall into four basic categories: halocarbon compounds, inert gas 
mixtures, water-mist or fogging systems, and powdered aerosols. The growing list of 
alternatives to Halon, in conjunction with advanced detection and fire resistant materials, 
provides protection from a broad spectrum of potential hazards. However, matching a 
fire-fighting agent to the specific uses required in confined areas on combatant ships poses 
a difficult problem. 

The Navy has attempted to identify non-ozone depleting alternatives to fill the role 
of HALON-1301 fire suppression systems. Many inert gas mixtures operate on the 
principal of oxygen displacement, quickly turning a confined space into a suffocation zone 
for escaping personnel. Dry chemical agents tend to be very corrosive to electronic 
equipment, causing significant damage in a shipboard environment. Traditional foam and 
water flooding systems are quite effective, but application must be limited due to stability 
considerations on a floating platform. Water mist systems seem to be the best-suited 
replacement for Halon flooding systems on U.S. Navy combatants. 

The U.S. Navy has developed a machinery space water mist system that utilizes a 
modified high-pressure spray nozzle. According to the Commission on Physical Sciences, 
Mathematics, and Applications (CPSMA) (1997), the system is designed to produce high 
volumes of 100 pm droplet (mean diameter) sprays with very high spray momentum to 
achieve rapid suppression of large hydrocarbon pool or spray (Class B) fires. These nozzles 
emit 2 gpm at 1000 psi on a grid system. This system produces water flow application rates 
of 0.06 — 0.07 gpm/ft 2 , three to four times higher than commercially available systems. 



Blast Mitigation using Water Mist 

The recent trend of using fine water mist systems to replace the legacy HALON- 

1301 fire suppression systems warrants further study into other applications of the water 

! 

mist systems. Preliminary research by the Naval Research Laboratories indicates that fine 
mists (20-25 pm droplet size) may reduce peak overpressures of a shock wave traveling 
through a space. Such pressure reductions could be used to mitigate the destructive effects 
of a shock wave (initiated by an explosive device) traveling through a structure. Although 
the mitigation effects are not observed to the same extent for droplet sizes in the 100 \xm 
range used in the current U.S. Navy water mist system, modification to smaller nozzle 
apertures such as those used on many commercially available water mist systems could 
make such a system feasible for installation on ships. 

Currently these blast mitigation effects have only been demonstrated in small-scale 
shock tube tests (Outa et al (1976) and Sommer field (1985)) and computer simulations 
(Schwer et al (2002 & 2003)). Uncertainty exists as to the scalability of such a system. The 
intention of this research was to fully investigate the applicability of such a blast mitigation 
system for shipboard use. A theoretical study of the mechanisms of blast mitigation using 
water mists was undertaken to explain the phenomena and possible methods of scaling and 
application. Current research into water mist blast mitigation was thoroughly reviewed and 
examined for trends to aid in guiding the direction of future work. In addition, a study of 
the degree of mitigation necessary to make a system practical for shipboard installation was 
conducted, including a close look at the calculation of failure loads for stiffened panels. 

These trends provide necessary information to begin design of a scalable blast 
mitigation system with applicability to real-world shipboard spaces. Given the recent trend 
toward tumblehome hull forms in future Naval Combatant designs, there exists strong 
applicability of this system in the “dead” spaces created by the shaping of the tumblehome 
hull. The end goal is a feasible design of a blast mitigation system to be used in the 
outermost spaces of Naval Combatants to protect interior vital system spaces. 






SHOCK STRUCTURES IN GAS-PARTICLE MIXTURES 
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The first phase of investigating the use of water mists for blast mitigation involves 
studying the dynamics of the mist-air mixture in the presence of a/blast or shock wave. In 
order to understand the phenomena observed in shock tube tests and numerical 
simulations, the governing dynamics of gas-particle mixtures were studied and applied to 
the specific case of a water mist. Although a complete theoretical explanation of the blast 
mitigation effects observed when a shock wave passes through a water mist - air mixture 
does not exist, the phenomena can be at least partially explained by further study into the 
dynamic interaction of pressure waves and gas-particle mixtures. 

The interaction of a shock wave with a mixture of water mist and air is a complex 
time and space dependent process. Understanding this interaction requires the use of 
simplifying assumptions to identify the properties of the gas mixture. The dynamics of gas- 
particle mixtures are heavily studied for low Mach number processes in works such as CKft 
et al (1978) and Crowe et al (1998). These studies employ the basics of fluid dynamics such 
as the Navier-Stokes equations for slow viscous flows, modified to include particle mixtures 
and interactions of heat and mass transfer for particles of various shapes and sizes. Crow et 
al (1998) develop particle interaction models, phase equations for multi-dimensional flows, 
and droplet-particle cloud equations from both a Lagrangian and an Eularian approach. 
These equations are then used as the background for numerical modeling techniques for a 
variety of multi-phase flows. 

Energy Methods 

Part of the numerical modeling involved in the above studies includes work to 
model evaporating droplets in sprinkler fire suppression systems. This includes energy 
equation modeling for individual droplets including convective and radiative heat transfer. 
Heat transfer could represent a significant portion of the total energy absorption 
encountered in the blast mitigation process, although the time dependency of the heat 
transfer process is likely to limit its effects when dealing with short duration blast pulses. It 
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is speculated that the time necessary to allow for heat transfer to occur to the mist particles 
is significantly greater than the time for the shock front to pass. The effects of heating and 
vaporization of the water mist are likely to play a major role in the cooling and fire 
suppression achieved, but not in the mitigation of blast overpressures. 

To further investigate the role of vaporization and heat transfer, an energy method 
was used to compare the energy required to heat and vaporize a unit volume of mist with 
the energy of the impinging shock wave. The heat of vaporization of fresh water is used to 
calculate representative energy absorption for vaporizing one cubic meter of water mist at 
room temperature. Appendix A contains the calculations demonstrating that for 25 jam 
water mist, using the U.S. Navy standard design mist concentration, 6 mega joules (MJ) of 
energy are required to vaporize one cubic meter of mist from room temperature. The 
energy of the shock wave is calculated by integrating the pressure of the wave over time and 
multiplying by the speed of wave propagation. In this case, the shock wave is assumed to 
propagate at the speed of sound. The energy of the representative shock wave used is 
0.239 MJ per square meter. 



Shock energy in 

Figure 2: Simplified Energy Flux for Unit Volume 

A net energy is calculated by multiplying the energy of the shock wave by the unit 

area over which it impinges and subtracting the energy of the unit volume of mist from this 

quantity. The result of this calculation is a negative energy, indicating that if the mist were 
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in fact vaporized to absorb the shock energy, the shock would be completely dissipated. 
Experimental results show that this is not the case, supporting the theory that the 
mitigation effects observed are not primarily caused by vaporization. Schwer et al (2003) 
provide further validation of the assumption that this vaporization energy is not a major 
contributor to the net mitigation effects through numerical simulations described later in 
this report. 

Gas-Particle Dynamics 

Investigation of the mechanics of shock propagation in an air-water mist mixture 
provides a background for understanding the mitigation achieved. Marble (1970) develops 
a relevant discussion of the mechanics of heterogeneous media for this study. In 
developing constitutive equations for mixed media, he makes a series of profound 
assumptions. Although each assumption is essential to the development of an analytical 
solution for a complex, mixed-media problem, it is likely that many of the phenomena 
causing the mitigation effects in water mist are eliminated from the analysis through these 
assumptions. Therefore, a careful look at Marble’s assumptions and their repercussions is 
essential to understanding his development and its applicability to the problem at hand. 

Marble’s first major assumption is that the volume fraction of the solid-particle or 
droplet cloud is considered to be so small that interaction between individual particles may 
be neglected or simplified. Since the mass loading experienced for fine water mist is 
relatively small, this assumption is not likely to invalidate Marble’s conclusions for this 
problem. Detailed transport processes between particles and gas are considered to be 
independent from the dynamic problem. Factors not considered in this analysis yet 
relevant to mist mitigation include the impact of particles on walls, and the effects of 
particle settling due to gravity. Marble considers the gas to behave as an ideal gas. 

A viscous stress tensor exists for gas only and is considered undisturbed by the 
presence of the particles. The volume fraction of the particles is considered negligible. 
This assumption allows omitting the Einstein correction from continuity equations. Marble 
uses a “smoothing 55 assumption for the continuum variables of the gas. Gas velocities and 





temperatures vary strongly in the immediate vicinity of a particle moving through and 
exchanging heat with the gas, but for the purposes of continuum equations, values are used 
which are averaged over a gas volume containing several particles. This assumption 
essentially mandates that all particle wakes and regions of immediate influence be dissipated 
very rapidly over the gas volume considered. The phenomena of particles forming groups 
and traveling in one another's wakes are not accounted for under these assumptions, 
although this grouping is not likely to play a major role in pressure reduction or energy 
absorption. 

To model the processes of drag and heat transfer for a coupled particle/gas cloud, 
the Reynolds and Mach numbers determine a regime for the flow and the corresponding 
governing equations. Generally, numerical solutions are required unless the flow falls into a 
clear-cut regime. such as Stokes flow, quadratic flow, or hypersonic flow. By limiting the 
analysis to a single regime, the limits of the approximation can be clearly defined and 
analytical solutions obtained. Marble assumes that Stokes flow and a Nusselt number of 
one describe particle drag and heat transfer. 

Appendix B shows the derivation of the governing equations for particles of radius 
a moving in a gas of viscosity p. Included in this derivation is a definition of the 
characteristic time, shown in Equation 1 below. This is the time required by a particle to 
reduce its velocity relative to the gas by e" 1 of its original value in an unaccelerated state. 

m 

Ty _ 

6 • 7i • a • |i (j) 

T v = characteristic time of particle 
m - mass of individual particle 
a = particle radius 
p = gas viscosity 

This equilibration time can be compared with the characteristic time (x) for the flow 
to further characterize the gas-particle interaction. For example, if x v »x, the particle 
moves through the flow field in a time scale short enough that the motion of the particle 
does not alter appreciably. In this case, the particle motion is dictated almost exclusively by 




its state before entering the system, independent of the properties of the gas. For the 
opposite case, when x v «X, the particle adapts to the local motion of the gas very quickly 
upon entering the system and moves only a very short distance before gaining equilibrium 
with the gas. In this case, the particle motion is relatively independent of its previous 
history and largely dependent on the motion of the local gas. 

The thermal behavior of the gas-particle mixture follows similar patterns to the 
kinetic phenomena described above. The structure of a normal shock wave exists in the 
region between complete equilibrium of gas and particles and the complete independence 
of gas and particles described above. Complete equilibrium exists both ahead of and 
behind the normal shock, a necessity of continuity. Therefore, conservation relations apply 
across the shock in terms of the modified gas properties for the gas-particle mixture 
described in Appendix B. The only deviation from this is that the existence of the shock 
structure is based on the upstream Mach number computed from the true gas properties, 
not the modified mixture properties. 

The shock structure itself consists of a gas-dynamic shock followed by a relatively 
thick zone in which equilibration processes take place. In order to separate these two 
zones, the relaxation time of the particles must be long compared to the molecular collision 
time of the gas. This assumption is reasonable since relating these two time quantities 
involves forming a ratio of the particle substance density to the density of the base gas. For 
fresh water mist in air, this ratio is high enough to allow significant leeway in the other 
factors involved such as particle radius and mean free path yet still meet the requirements 
of the assumption. 

Marble postulates that the gas-dynamic shock structure is unaffected by the particle 
cloud, though the shock may be thin relative to the particle radius. Because the volume 
fraction of the particles is negligible, at any time less than 1% of the shock front is being 
punctured by particles. Therefore, the particles pass through the shock front unaffected by 
the shock and retain their initial conditions upon entering the equilibration zone. 





It is in this equilibration zone where Marble postulates the particles have their 
greatest effect. The relaxation lengths (obtained by integrating particle speed over the 
relaxation time) for both velocity and temperature determine the thickness of the 
equilibration zone. The thickness of this zone subsequently determines the shape of the 
exponential decay following the shock front. The generalized shape behaves as exp[- 
x/(M 3 A)] where M 3 is the localized Mach number in the equilibration zone and A is the 

total relaxation length of the gas-particle mixture for either velocity or temperature 

) 

equilibration. Solutions for A include a first order dependence on the relaxation length of 
the particles and a series of terms including a factor of r|, the particle mass loading, which 
represent the gas-particle interaction. Figure 3 shows the spatial structure of the relaxation 
zone. It can be observed that the effect of the particles essentially extends the equilibration 
zone, allowing more time and space to dissipate the shock energy. 



Figure 3: Spatial Structure of Shock Relaxation Zone (Marble (1970)) 

To this point, no consideration is given to the effect of different sized particles and 
collisions between particles. As long as the expected time between collisions of particles is 
equal to or longer than the equilibration time for the particles, a particle will essentially lose 
the scattering velocity obtained through a collision encounter before it experiences another 
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collision. This assumption simplifies the problem by always considering the starting 
conditions to be equivalent to the equilibrium slip condition. 

Sommerfield (1985) discusses the decay of incident shock waves upon encountering 
gas-particle mixtures from both a velocity and a peak pressure standpoint. The general 
premise is that shock waves initiated within or propagating into a gas-particle mixture will 
have differing characteristics from the pure gas case due to the inertia and heat capacity of 
the particles. There exists a transition region as a shock wave encounters a gas-particle 
mixture where the velocity of the shock wave changes from the initial value to a lower, 
equilibrium value. 

Building on the work of Marble (1970), Sommerfield also develops “equivalent gas 55 
properties for the gas-particle mixture based on the loading ratio (by mass and specific heat) 
of particles to gas. Depending on the properties of the gas-particle mixture, the shock wave 
is refracted and reflected (as a shock or a rarefaction wave) at the interface. For the case of 
ship geometry and a water mist filled space, this interface coincides with a breached 
bulkhead or portion of ship hull, so the shock front would also be modified due to the 
geometry of the hull as it deforms. 

The velocity equilibration zone following a shock is extended by the presence of a 
gas-particle mixture. Due to the heat capacity and inertia of the particle or droplets in the 
mixed media, the velocity of a transmitted shock wave does not change instantaneously. A 
transition region develops where the shock velocity is slowly reduced from its initial to its 
equilibrium value obtained from the equivalent gas analysis. During this deceleration 
process, peak pressures in the shock front are increased slightly due to basic Bernoulli flow 
considerations. As velocity decreases, pressure must increase to maintain equilibrium. This 
pressure increase occurs before equilibration of the gas-particle mixture begins to dissipate 
the energy of the shock front. 

In terms of mitigation effects, the reduction in velocity represents a reduction in the 
overall strength or energy of the shock front. Although it is not instantaneous, if the 



velocity of the shock is reduced prior to the next impingement on a structure, the energy of 
the wave on this structure is effectively reduced. However, enough time and space must 
elapse to allow the enhancement of peak pressures due to the shock deceleration to pass 
before mitigation effects take over and peak pressures are reduced. The depth of mist 
necessary to achieve mitigation and not enhancement must be quantitatively defined to 
properly design a mitigation system. 

These theoretical developments help to explain the phenomena recorded from 
Outa et al (1976) and Sommerfield (1985) from shock tube testing. The results of these 
tests and numerical simulations by Schwer et al (2002 & 2003) provide background of the 
mitigative effects of water mist. The next section studies previously published experimental 
and numerical results and places them in the context of the water mist-air mixture behavior. 

Experimental Results 

Outa et al (1976) measure the time history of a shock wave initiated in the pure gas 
section of a shock tube and subsequently interacting with a dusty gas and decaying through 
the mixture until equilibrium is achieved. A graphical representation of the decay of the 
shock wave as it moves through the heterogeneous zone of a shock tube test is shown in 
Figure 4. The bend in the shock front on this time versus distance plot represents the 
slowing of the wave propagation as the mixed media is encountered. 



Figure 4: Wave Diagram of Shock Tube Flow Showing Decay of Incident 
Shock Front Through Gas-Particle Mixture (Outa et al (1976)) 
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A particle loading ratio of less than two consisting of 3-9 pm diameter white carbon 
particles are used in these shock tube tests. The loading ratio for the current U.S. Navy 
water mist fire suppression system is 2.0. Therefore, this loading ratio has good 
applicability to the water mist problem, although the mist droplets are significantly larger 
than these particles. Through testing and theoretical calculations, Outa establishes trends 
for the effects of particle diameter on pressure, temperature and velocity. These results are 
shown in Figure 5 below. The pressure plot in this figure is the most relevant to this 
problem. Particle sizes below 15 pm appear to delay (in a spatial sense) the arrival of the 
peak pressure but do not appear to cause a significant decrease in the peak pressure itself. 
However, for the 30 pm particle size, the pressure equilibrium value seems to be lower than 
for the other cases, indicating a mitigation effect. 



Figure 5: Effect of Particle Diameter on Temperature, Velocity, and Pressure 
(Outa et al (1976)) 

Outa conducted shock tube tests for various mass loadings. Figure 6 shows plots 

of the shock as it passes various locations along the length of the seeded section of the 

shock tube. The traces show a discontinuity in the leading edge of the pressure trace as the 

shock front passes the sensor. This is representative of the steep pressure increase at the 

start of the shock wave. Comparing the traces along the 2 meters of shock tube accounted 

for by the various traces shows the weakened pressure wave passing through the gas- 

particle mixture. The mitigating effects of the particles are obvious from the pressure drops 
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observed between the four traces in each plot series. The total mitigation due to the mixed 
media is represented by the difference in the pressure between the upper left-most trace 
and the lower right-most trace in each row of plots. These traces are from the x locations 
corresponding to the start and end of the mixed-gas section of the shock tube. The three 
rows represent three different mass loading cases at slightly different Mach numbers. The 
low mass flow loading case (r] 0 = 0.6, top plot) shows the least mitigation, reducing the 
peak pressure by about 20% over the length of the test section. In contrast, the T] 0 = 2.1 
case shows a reduction of the unmitigated shock pressure by 30%. 
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Figure 6: Pressure Traces Showing Variations of Wave Form Along Tube 
(Outa et al (1976)) 


To better demonstrate the pressure reductions achieved by the gas-particle mixture, 
Outa plots the experimental results of the shock pressure against the expected, unmitigated 
shock structure. Figure 7 shows the experimental pressure trace and calculated stationary 
wave forms at a section where the leading front is still decaying. The experimental 
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pressure, as indicated by the open circles, approaches an equilibrium value behind the 
shock front. This pressure is significantly lower than those calculated for stationary waves. 
The peak pressure is again reduced by about 30% from the expected pressure prior to the 
shock encountering the gas-particle mixture. 



Figure 7: Experimental Pressure Trace and Calculated Wave Forms (Outa et al (1976)) 

Sommerfield (1985) repeats many of the experimental results presented above using 
a similar shock tube arrangement These experiments were conducted using glass spheres 
with a median particle diameter of 27 pm rather than carbon dust as above. The particle 
size used by Sommerfield duplicates a feasible droplet size for water mist systems using 
current nozzle technology. The droplet size for a mitigation system must be optimized to 
achieve the desired mitigation effects for a given system design. The choice of droplet size 
for the proposed mitigation system will be discussed later in this report. However, it is 
noted at this point that the particle size used in Sommerfield’s experiments corresponds 
directly to the proposed droplet size for a water mist mitigation system. In addition, the 
density of the glass spheres used is not significantly different from the density of a fresh 
water mist droplet. The resulting pressure traces, shown in Figures 8 and 9 below, clearly 
demonstrate the mitigating effects of the particle mixture used by Sommerfield. 














Figure 8: Pressure Traces for M so = 1.49, rj = 0.63 (Sommerfield (1985)) 

The first experimental results shown are for a mass loading ratio of 0.63. Figure 8 
shows a total of 8 pressure traces, representing 8 locations along the seeded portion of the 
shock tube. The numbers on the left side of the trace show the zero mark (pre-shock 
pressure) for each location as a reference point. The pressure jump as the shock passes 
causes a sharp discontinuity in the pressure trace from this initial value to the more readily 
observed trace above. The first notable feature on the plot is the decay of the sharp shock 
front as it moves through the gas-particle mixture. This decay is even more evident in 
Figure 9, the results of a similar experiment at a mass loading ratio of 1.33, over double the 
loading shown in the previous result. 
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Figure 9: Pressure Traces for M s0 = 1.127, r| = 1.33 (Sommerfield (1985)) 

The decay of the shock front observed in these traces represents the effects of the 
deceleration of the shock wave as it encounters the gas-particle mixture. The net effect of 
this velocity reduction is observed in the pressure traces as a delay in the arrival of the peak 
overpressure from the shock. Rather than an instantaneous increase to peak pressure 
followed by decay as observed in a pure gas, the pressure increases more gradually in the 
seeded gas, resulting in a delay in the arrival of the peak pressure. 

The distance between the first and last pressure gauges recorded in these plots is 
approximately 4 meters. The peak pressure between these traces is reduced by 28.6% for 
the lower mass loading case. A reduction by 20% results between the first and the 4 th 
sensor, a distance of 1.25 meters. For the high mass loading case, the total reduction over 
the length of the tube is slightly lower than that for the lower mass loading case. The decay 
of the shape of the shock front is greater for this case, but it appears that the lower mass 
loading provides better overall pressure mitigation characteristics. This observation 
correlates to the theoretical development of the deceleration zone presented in the previous 
section. The higher mass loading causes more deceleration of the shock front, extending 
this zone of the overall structure and therefore extending the pressure enhancing zone and 
leaving less time for pressure mitigation to occur before the system reaches equilibrium. 
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Compared to the ratios studied by Sommerfield (rj—0.63-1.33), current IJ.S. Navy 
water mist systems provide higher mass loadings (r|=2.0) of water. Scaling this mass 
loading for blast mitigation purposes without sacrificing the back fit capability to current 
Navy systems is an important factor in making such a system feasible. Recognizing that the 
mass loading for the current Navy system is much greater than most commercial systems 
using smaller droplet sizes, the mass loading is recalculated based on a smaller 
concentration of water mist. The concentration is reduced to Vi of the current Navy design 
based on data from commercial water mist systems and the reduction of the mean mist 
diameter applied by a factor of four from the baseline Navy design. This concentration is 
more representative of the mist applied through the smaller aperture nozzles needed for 
this droplet size. For the purposes of calculating the behavior of shock waves through a 
representative mist-filled space, a mass loading ratio of 0.512 is used. 

Based on these experimental and theoretical calculations, mitigation on the order of 
30% reduction of the peak overpressure of an impinging shock should be possible. Proper 
droplet size and mass loading of the water mist are essential to achieving this mitigation. 
For further investigation into the mitigation possibilities of water mists, numerical 
simulations must be used. 




MITIGATION EFFECTS 


As observed by Schwer and Kailasanath (2003), shock tube tests leave many 
unanswered questions about the behavior of seeded gases and shock fronts. In 
experiments where the driven section and the driver section of the shock tube are of similar 
size, the mitigation effects are small. The primary effects observed are a reduction in the 
velocity of the shock front, and a reduction of the pressure at the shock front. The overall 
maximum overpressure occurs slightly behind the shock front due to the deceleration of 
the shock and the equilibration zone behind the shock. The maximum overpressure, 
occurring further back in the shock structure, is initially increased in computer simulations 
of this setup. This effect is due to the “compression” of the shock front due to its 
deceleration upon impingement with the seeded gas. 

The constitutive equations presented by Outa et al (1976) and Sommerfield (1985) 
were adapted to the specific case of the air-water mist system by Schwer et al (2003). The 
result of these calculations is an accurate numerical modeling tool for use in predicting the 
mitigation effects and behavior of a shock wave passing through an air-water mist mixture. 
Integrated into this simulation are the relevant factors to simulate the use of liquid droplets 
rather than the dry particles used in all laboratory testing to date. 

Results from Schwer’s simulations indicate that overall mitigation is tied to the 
location of the maximum overpressure. For shock waves traveling through seeded gases, 
the maximum overpressure location no longer coincides with the shock front location due 
to deceleration effects. As shown in the experimental results above, at higher mass loading 
ratios the sharp shock front decays, resulting in a delay in the arrival of the peak 
overpressure. From Schwer’s numerical model. Figure 10 provides a distance-time plot 
showing the arrival of the shock front and the arrival of the peak pressure for a variety of 
particle sizes at a mass loading of 0.5. The split between location of the shock front and 
location of the peak overpressure is closely correlated to the particle size. Larger particles 
cause a longer delay of the arrival of the peak pressure in these numerical simulation results. 



Schwer’s numerical simulation results shown are well correlated to the pressure traces 
shown above. The differential in arrival time of the peak overpressure is easily observed in 
Figure 9, although the correlation presented by the data in Figures 8 and 9 relate to mass 
loading ratio, not to particle size. These two factors exhibit the largest influence on the 
mitigation characteristics of a mixed-media and therefore are the main design variables for 
creating a feasible mitigation system using water mist. 



Figure 10: Location of Shock Front (solid) and Maximum Overpressure 
(dashed) for Different Particle Sizes (Schwer et al (2003)) 

Although the effects of the seeded gas on the shock wave velocity and shock front 
pressure rise are obvious from these experimental and numerical results, the maximum 
overpressure initially increases for many cases simulated using this shock tube arrangement. 
It is theorized by Schwer et al (2003) that the relative sizes of the driver and driven sections 
of the shock tube arrangement are the root cause of this increase. The setup used has a 
driver section 3.5 meters in length and a driven section 4.0 meters in length, simulating a 
blast front of the same size as the medium through which it travels. However, explosive 
shocks act much more like a point load in a large medium. The driver section of a real- 
world explosion is significantly smaller than the surrounding, “driven” section, which 
experiences the effects of the shock front. Therefore, a more realistic model was simulated 
numerically by Schwer using a driver section 1/20* the size of the driven section. 






The simulations run by Schwer with this more realistic geometry demonstrate a 
closer correlation to the actual mitigation levels achieved in a laboratory setting. The 
maximum overpressure was measured for a variety of particle sizes and plotted as a 
function of distance along the seeded gas portion of the tube. The effects of droplet size 
on peak overpressure for a constant mass loading ratio are shown in Figure 11. From the 
plot on die left, the optimal droplet size appears to be between 20 and 40 pm, with the 40 
pm droplet size beginning to exhibit a lag in the equilibration time behind the shock front. 
The plot on the right demonstrates that the mitigating effects of larger droplets (100 pm) 
such as those used in the current Navy design are not as pronounced as for the smaller 
droplet size. 



Figure 11: Snapshots of Pressure Profiles as 5 ms intervals for M D = 1.49, r\ = 

0.63 (Schwer et al (2002)) 

The peak pressure observed just as the shock front crosses into the mixed gas area 
is slighdy higher than the no-particles case due to the deceleration of the shock front. 
However, the reduction in peak pressure as the shock front travels through the seeded gas 
is dramatic. For most particle sizes, the majority of the mitigation of pressure has occurred 
after the shock has traveled through 5 meters of the seeded mixture. The level of 
mitigation itself is dependent on particle size. Mitigation appears to increase with particle 
size to a maximum for the 25 pm particles. The results for the 60 pm particles show a 








slower decay to the maximum mitigation level, although the pressure at a distance of 8 
meters is nearly equal to that achieved for the 25 pm particles. 

In comparing the results presented in Figure 10 to those presented in Figure 11, it 
appears that the overall mitigation is tied to the location of the maximum overpressure. 
The mitigation is achieved by a reflected rarefraction wave dissipating the shock pressure. 
Situations where the maximum overpressure is farther behind the shock front result in the 
reflected rarefraction wave dissipating the pressure sooner, allowing for greater overall 
mitigation. To some extent, larger droplet sizes achieve greater mitigation under these 
assumptions. However, these results are for a limited, shock tube geometry. The 
optimizations are still based on geometry and mass loading of the real-world space through 
which the shock wave travels and the shape and size of the initial shock. In order to 
investigate the dependencies of these factors, a more realistic, spherical blast equivalent to 
4.67 lb of TNT was simulated by Schwer et al (2003). 

The results of Schwer’s simulation of this geometry for unmitigated and mitigated 
blasts on a time basis are shown in Figure 12. These simulations were conducted for the 
same mass loading (0.5) as above and for an initial droplet size of 25-30 pm. Both the time 
delay and pressure reduction effects of the water mist impingement are observed. The 
mitigation effect is not observed at the earliest time interval (125 psec), but is obvious by 
the next increment (375 psec). This plot also shows the travel of the shock wave 
corresponding to the time intervals. The first mitigation effects are observed on the second 
time interval, at a distance of 0.8 meter from the initial encounter with the seeded gas. 
Although this model represents a geometry not validated by laboratory experiments, it does 
provide a rough estimate of the depth of mist necessary to achieve mitigation effects. 
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Figure 12: Blast Pressure Profile at Different Times for Unmitigated (Solid) 
and Mitigated (Dashed) Tests (Schwer et al (2003)) 

Vaporization effects were also considered as part of Schwer’s numerical 
simulations. The results of the simulations support the theoretical calculations from 
Appendix A, namely that vaporization plays a minimal role in the overall mitigation effects 
of the water mist. Figure 13 shows Schwer’s comparison of the peak overpressures for the 
pure gas, the seeded gas without vaporization effects, and the seeded gas including 
vaporization effects. The vaporization effect does not appear at all until the shock has 
traveled about a meter through the seeded gas. This corresponds to the minimum time 
required for any heat transfer to occur from the blast wave. Although this distance also 
corresponds to the depth of mist required to overcome the pressure enhancement zone due 
to shock deceleration, the relatively small deviation from vaporization shown in Figure 13 
supports the theory that the effects observed are not due primarily to vaporization. The 
overall blast pressure is reduced somewhat due to vaporization, but the majority of the 
reduction is due to momentum and energy transfer and not heat transfer and vaporization. 
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BLAST RESISTANCE OF STIFFENED PANELS 


Blast effects are intimately connected to the structural design of warships. Navy 
design guidelines lay forth requirements to withstand certain levels of static and dynamic 
loading on various portions of a ship structure. Although resistance to blast from outside 
sources other than underwater and nuclear explosion is not strictly defined for all U.S. 
Navy ships, a fundamental and detailed understanding of shock wave loading and the 
associated structural response of the ship is required for effective warship design. In 
addition to aiding in proper design of warships, a full understanding of the behavior of ship 
structures to blast loading defines the role that a mitigation system might play. The level of 
mitigation necessary to make such a system practical is determined by the reduction in 
loading required to maintain the integrity of the ship. 

Analysis of the behavior of a floating, multi-degree-of-freedom system such as a ship 
under blast loading is best accomplished through study of live fire data. The late 1940 s 
produced a number of landmark studies into the failure of ships under explosive force. 
The experiences of two World Wars provided ample data on torpedo attack and 
underwater explosions. Taylor (1948) and Cole (1948) conducted comprehensive analyses 
of submerged blast waves and their effects on thin plates. Their research formed the 
technical basis of all current blast damage prediction methods. Although most subsequent 
research into this failure mode, conducted by ONR, has been classified confidential, open- 
source study has been conducted on blast loading of steel structures and stiffened panels. 

Nonaka (2000) looked at shear failure of steel beams and concrete structures using the 
World Trade Center bombing in February 1993 as an example. Rudrapatna et al (1999 and 
2000) conducted blast loading tests of square plates and small, stiffened plates. These tests 
provided background for the possible deflection modes of a larger scale stiffened panels. 
Shin et al (2003) defined design guidelines for explosive confinements using analysis of 
shock transmission and reflection, important to understanding the proper loading to apply 
in the context of this problem. Schleyer et al (2000) provided background on numerical 



modeling of pulse pressure loading. Ramajeyathilagam et al (2000) provided shock pressure 
pulse loading models including equations for the shape of the time dependent pressure 
pulse, peak pressure, and decay constant based on the weight of explosive. Nurick et al 
(1995) studied deformation and tearing for small square plates with a single stiffener milled 
into the test sample. Again, this study provided a good idea of mode shapes for larger scale 
stiffened panel deflection. Klaus (1985) looked at numerical modeling of panel walls 
subject to blast loading. Finally, Houlston et al (1985) developed and validated finite 
element models for blast loading of stiffened panels by comparison to experimental work 
done on blast loaded square plates. Houlston’s work is quite relevant to this study, 
although only carries the analysis to plastic deflections, not to full failure of a panel. 

The behaviors of plates, stiffened panels, and other steel structures have been 
extensively documented through this substantial body of literature. Ship structures are 
comprised of stiffened panels, composed of steel plate with evenly spaced steel stiffeners in 
the longitudinal direction. Larger structural members such as girders and frames support 
each panel. For the particular case looked at here, the confining end conditions will be at 
the frames and girders and the panel affected will be a stiffened panel with uniform size 
stiffeners. These assumptions are reasonable considering the fabrication of typical naval 
combatants. Specific panel geometries installed on representative ships are presented in a 
later section of this report. 

When transverse loading of any kind is applied to such a stiffened panel, the panel 
will deflect to support the load as expected by statics and Hooke’s Law. The dynamics of 
blast loading conditions modify the behavior of the panel; time dependencies and impulse 
effects become factors in the deformation and yield characteristics of the panel. Regardless 
of the dynamics of the problem, the yield behavior of the material will have certain 
ramifications on the panel deformation. 

As part of the panel yields, it can no longer support shear stresses and bending 
moments. The stresses and moments on the panel are therefore transferred outward to the 
non-yielding portion of the panel until the energy of the blast wave is absorbed or reflected. 




The yield zone continues to grow until this equilibrium is achieved. Stiffeners or girders 
will slow the spread of the yielding zone since they will require greater stress concentrations 
to yield and will absorb more strain energy in yielding. 

Despite understanding the yield mechanisms and dependencies of the problem, no 
simple, reliable method of predicting hull plate blast damage has been developed. Although 
computer codes for the prediction of blast damage are available, most codes provide 
prediction of elastic and plastic deformation only. Prediction of fracture onset and 
propagation is not well validated against actual results, so the codes provide little more than 
a rough estimate of potential damage. As a result, nearly all of the definitive blast damage 
prediction is conducted using scaled, live fire tests, requiring substantial time and resources. 
For the purposes of this study, prediction of the overpressure required to fail a bulkhead or 
boundary is necessary. In order to make such a prediction, “failure” itself must be 
quantitatively defined. 

i 

Once the definition of failure is understood, the loading that causes this level of 
damage must be defined. Prediction of the level of blast loading a structure can withstand 
is central to understanding the level of mitigation necessary to preserve the integrity of the 
structure. To establish a clear criterion for the level of mitigation to be achieved, a 
representative blast load must be analyzed and compared to the loading that the ship 
structure can resist. The shock loading created by various explosives is well defined in 
open-source literature. Therefore, with information on loading and failure mechanisms of 
ship structures and data on the potential blast loads encountered, the level of mitigation 
that must occur to prevent failure can be defined. 

A confined space mitigation system such as water mist is not likely to prevent 
damage to the outer hull of a ship. The goal is to protect the next layer of watertight 
boundaries in the ship to prevent the spread of damage. Confined space geometries 
provide added challenges to determining the actual level of blast overpressure propagating 
through the space. This phenomenon has been studied in open-source literature for blast 
mitigation in civilian, land-based structures. Explosions confined by a chamber or room 






result in rapidly increasing gas pressures, which eventually decay by venting out of either a 
newly breached bulkhead or through the entry point. Under these conditions, shock 
reflections occur such that the overall effect can be greater than the incident shock. 

A study must be undertaken as to the dynamic pressure (blast) loading that typical 
U.S. Navy ship’s interior bulkheads and decks can withstand. The underlying assumption 
in this study is that preserving the integrity of a bulkhead or deck minimizes the spread of 
damage past that boundary. Prevention of first yield, however, is a very conservative, 
unrealistic expectation when dealing with blast loading. Mere plastic deformation does not 
define failure in this sense. Significant (macroscopic) tearing or cracking of the plate or 
weldment material will be considered failure for the purposes of this study. Under this 
definition, the boundary in question will no longer be watertight or structurally sound. Full 
buckling of a panel would also be considered failure since the boundary could no longer 
provide resistance to normal ship loads. Buckling in this sense is likely accompanied by 
cracks or fractures along the panel edges and therefore is included in the definition above. 

Blast Load Prediction 

The first step in determining the level of mitigation necessary to preserve the 
integrity of boundaries is to determine the blast profile to be resisted. This profile consists 
of three major components: duration of overpressure phase, peak overpressure, and 
impulse magnitude. For explosions close to the target, pressure-driven effects occur 
quickly. The characteristic duration of a high explosive detonation is measured in 
microseconds (10' 6 ). The pressure-driven effects of such an explosion occur on similar 
time scales. 

Air-blast loading associated with explosions is commonly subdivided into (1) 
loading due to the impinging shock front, its reflections, and the greatly increased static 
pressure behind the front, all commonly known as overpressure; and (2) the dynamic 
pressures due to the particle velocity and mass transfer of the air. Pressure loadings are 
customarily characterized in terms of scaled range, given by: 






( 2 ) 


Z = R/W 1/3 
Z = scaled range 

R = radial distance between explosion center and target 
W = explosive weight (TNT equivalent) 

Units for charge weight and distance in this expression should be either pounds and feet or 

kilograms and meters. Using this scaled range factor, constant values of Z result in 

constant parameters for explosive effects (peak pressure, positive duration, etc.). 

The explosive effects parameters such as peak pressure and positive duration are 
dependent on the size and type of the explosive. However, characterization of the shape of 
the blast profile is somewhat standard and often parameterized. For the purposes of this 
study, the preliminary estimates of the blast loading encountered by USS COLE will be 
used as a representative case. 

Appendix C contains calculations for a 500 kilogram TNT equivalent explosion at a 
standoff distance of 10 meters using various methods for parameterizing the blast loading. 
The time and distance dependent shape of the actual pressure profile of the shock wave is 
an important parameter to understanding the blast loading to apply to a panel. The shape 
of the incident wave is given by Taylor (1941) as: 

Pi = P G e fl(t - x/c) (3) 

Pi = pressure in incident wave 
Po = peak overpressure 
t = time elapsed 

x = distance from blast center, measured perpendicular to impinging plate 
This shape reflects an exponential decay from the peak overpressure based on both time 
elapsed since the blast and distance from the blast center. The empirical constant n in 
Equation 3 is a decay constant representing the nature of the decay. Ramajeyathilagam et al 
(2000) give an empirical equation for calculating the inverse of this decay constant (© = 
1/n) as follows: 
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W = charge weight (kg TNT) 

R — radial stand-off distance (m) 

The peak overpressure used to determine the shape of the incident wave can be 
calculated using empirical formulae or charts as shown in Appendix C. Nonaka (2000) 
calculates peak overpressure using a log-based empirical formula and an equivalent mass of 
Sinkiri Dynamite. The resulting overpressure for the representative 500 lb TNT explosive 
is 130.3 psi or 89,850 Pa. Remennikov, A. (2004) uses the scaled range (see Equation 2, 
above) as an input to a graphical representation of various blast parameters (Figure 14). 
The resulting overpressure is approximately 145 psi or 10 6 Pa. Given the error inherent in 
reading values off the log-log graph shown in Figure 14, these results are remarkably close 
to those calculated using Nonaka’s empirical formula. Ramajeyathilagam, K. et al (2000) 
use an empirical constant multiplied by the cube root of the scaled range (Z). This 
calculation returns an overpressure of 7000 psi, over two orders of magnitude greater than 
Nonaka’s formulation. It is believed that there is an error in the empirical constant used in 
this formulation, resulting in this drastic difference. 
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Figure 14: Blast Parameters Based on Scaled Distance 
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For the purposes of this report, the overpressure used to calculate loading on ship 
structures will be 130.3 psi with a positive phase pulse duration of 0.2 thousandths of a 
second. To simplify the loading used for theoretical calculations, this pulse was assumed to 
be a square pulse of constant pressure over the positive phase duration. Although this 
square pulse is an accurate model of the leading edge of a normal shock structure, the decay 
experienced in this sharp discontinuity due to the mist mixture is not accounted for. The 
effect of the shape change in the shock front is minimal on the way the shock is applied to 
the structure for analysis. The dynamic effects modeled in this study accounted for a short 
duration pulse at peak pressure. The time duration of the pulse itself is not significantly 
altered by the change in shape of the shock, so the effects of this change are negligible. 

The shock described is the loading desired to be resisted by a ship structure for the 
purposes of this report. The remaining problem in determining the mitigation level needed 
is to define what level of shock will cause failure of a ship’s structure. The first step in 
developing an answer to this question is to understand how ship subdivisions are designed 
in current practice. 

Naval Design Standards 

When the structural arrangements of ships are designed, bulkheads and decks are 

designed to withstand the typical global and local loads applied to them without yield. 

These loads are due to large-scale ship motions such as hogging and sagging, weather 

conditions such as icing and green water, hydrostatic pressure on the hull, and the weight of 

equipment and payloads. Combatants must also comply with a series of Design Data 

Sheets (DDS’s) and with the General Specifications Tor Ships of the United States Navy 

(Gen Specs). DDS’s are design guidelines that outline a systems-based approach to the 

mitigation of damage. Specific DDS’s relevant to ship survivability and structural design 

include: DDS 072-1 (1972) “Shock design values/ 5 DDS 072-3 (1988) “Conventional 

weapons protection (fragments)/’ DDS 072-4 (1986) “Hull, mechanical, and electrical 

systems survivability,” DDS 072-6 (1987) “Shaped charge warhead weapon effects data,” 

DDS 072-7 (1988) “Conventional airblast (proximity),” DDS 072-8 (1986) “Conventional 

airblast (contact and internal) design and analysis methodology,” DDS 079-1 (1975) 
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“Stability and buoyancy of U.S. Naval surface ships/ 5 DDS 100-4 (1982) “Strength of 
structural members/ 5 DDS 100-6 (1987) “Longitudinal strength calculation/ 5 DDS 100-7 
(1984) “Structure to resist weapons firing effects/ 5 and DDS 100-9 (1991) “Nuclear airblast 
design for surface ship structures. 55 DDS 5 s 072-1, 072-3, 072-6, 072-7, and 072-8 are 
classified Confidential and therefore cannot be further discussed in this study. They set 
forth some of the major design criteria for blast loading on ship structures and therefore in 
many respects define the strength of ship boundaries. Although the design guidelines 
themselves are classified, the resulting structures of current surface combatants provide a 
good example of the type of structure that fulfills the requirements of the DDS 5 s. 

Most of the requirements set forth in the unclassified design guidelines are general 
in nature. For example, DDS 072-4 defines three major design principles for HM&E 
survivability: separation of redundant systems, isolation of non-redundant systems, and 
concentration, localization and shielding/armoring of vital non-redundant systems from 
areas of potential hazard. It then lays forth a 5-step procedure for designing a total system 
for redundancy and separation, two of the key survivability characteristics. DDS 079-1 
contains key parameters for proper stability characteristics dependent on ship type and 
mission. It defines requirements for parameters such as intact and damaged heel and trim 
angles and righting arm curves in various conditions. 

DDS 100-4 characterizes some important factors for stiffened panel calculations 
that are relevant to this study. First, it delineates a method of calculating buckling loads of 
plating panels under combined compressive and shear loading. It also identifies the Navy 
standard strength curves for selected common hull materials. These curves include column 
strength, ultimate compressive strength, and strength under combined compressive and 
shear loads for high strength steel (HS) and ordinary steel (OS). 

DDS 100-7 and the Gen Specs define requirements for hull structures based on 
the proximity of the structure to own ship’s weapons systems. For example, these 
guidelines set forth an empirical formula relating a static equivalent pressure to the booster 
motor thrust, angle of incidence, and area of impingement for missile blast loading from a 
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missile being launched from own-ship launchers. The intention of these guidelines was to 
protect ship bulkheads from the blast loads experienced behind the Mark 26 deck-mounted 
missile launcher. Although these launchers have been phased out of US Naval platforms in 
favor of the in-hull Vertical Launch System (VLS), the empirical calculation of equivalent 
static head pressure based on rocket blast loading is a useful tool. The DDS defines the 
static equivalent pressure applied by a missile launch as: 

T*[ sin(a) + 

^ P = —^- 

A (5) 

P = static equivalent pressure (psi) 

T = total thrust of missile (lbs) 
a= angle of incidence (deg) 

A = impingement area of surface (in 2 ) 

Equation 5 is not a direct function of distance from the missile to the bulkhead or 
deck in question; the area of impingement is calculated based on the projected area of a 3° 
divergence blast cone from the axis of the missile to the surface. Therefore, the distance is 
accounted for by documenting the impingement area and angle of incidence of the blast 
load to the surface in question. As the distance from the missile increases, this projected 
area also increases, reducing the static equivalent pressure experienced. 

DDS 100-7 also provides guidelines for plate panel design to withstand the 
equivalent pressures of missile and gun blasts. Empirical formulas and tabulated values of 
the relevant constants are provided for two loading cases. The formulas are developed 
assuming simply supported edges on the plates. Plate dimensions are determined by the 
locations of stiffeners and girders. The first loading case assumes that the blast load is 
uniform over the entire plate, indicating that the impingement area is larger than the area of 
plating between stiffeners. The second loading case assumes a more concentrated load in a 
central rectangular area of the larger plate section. 

Similar expressions are developed for gun blast loads based on the geometry from 
the structure being assessed to the gun muzzle and on the caliber of the gun itself. DDS 






100-7 also delineates a method for design of deck and bulkhead stiffeners once the loading 
is understood. It treats the stiffeners as pinned-pinned beams, therefore rendering panel 
stiffeners as non-continuous members. If panel stiffeners are considered continuous, some 
level of fixity could be assumed for the end conditions. The maximum moment in the 
stiffener is calculated assuming either a uniform or a concentrated load condition. Under a 
uniform load condition, the required section modulus equals the maximum moment over 
the maximum allowable stress. This maximum moment is calculated as shown in 
Equations 6 and 7. 

W = P n ~ (6) 



8 (7) 

W = distributed load on stiffener (Lbf) 

P = static equivalent blast pressure (psi) 
n = length of minor axis of impingement area (in) 

L = stiffener length (in) 

For a concentrated load case, the force applied is calculated as the static equivalent 
pressure times the area over which it is applied, similar to the formulation in Equation 6. 
This force is then applied as a point force at a pre-determined location along the stiffener. 
The maximum moment is defined using the distances from the load to each of the pinned 
ends times the force over the length. This point load method is a very conservative 
approach to design where the blast impingement area is significantly less than the stiffener 
length. It is much more common to assume the blast to be uniformly distributed over the 
stiffener for blast loadings. 

Design of major supporting structures such as girders follow similar procedures to 
those for stiffeners. Again, structural members are conservatively assumed to have pinned 
connections at support points and are subjected to either uniform or concentrated loads 
depending on the geometry of the blast zone. 

DDS 100-9 defines a procedure for “the design of surface ship structure to resist 
in-air, non-contact tactical nuclear weapons explosions of low- to medium-yield. 55 This 
design guideline is based on a blast and shock wave originating just above the water surface 
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with the blast wave front perpendicular to the water surface. The ship structure, including 
intermediate stiffeners between girders, frames, and transverse bulkheads, is designed to 
support the reflected overpressure from a blast and to resist the shearing, tearing, buckling, 
or collapse loads that exist on various faces of a structure as a blast load passes over it. The 
reflected overpressure is the resulting pressure occurring when the incident overpressure 
wave strikes a surface parallel to the blast front and is reflected. This pressure includes the 
incident overpressure plus the instantaneous pressure loading of the reflected shock front. 
DDS 100-9 uses a closed box as its global structural model and therefore looks at loads 
resulting on all outward faces of this structure as the blast travels over and around it. The 
maximum loading will occur when a blast impinges perpendicular to a panel. This loading 
represents the relevant load case for the stiffened panel problem presented in this’chap ter. 
Although this report primarily focuses on conventional blast loading, the shock wave 
formation of a conventional blast is not significantly different from that of a larger, nuclear 
blast at greater stand-off distances. 

Empirical formulas are provided in DDS 100-9 for the peak reflected pressure and 
stagnation pressure a structure must withstand. These pressures are calculated based on 
the peak incident overpressure in the shock and the ambient pressure of the surrounding 
atmosphere as shown in Equations 8a, 8b, and 8c. 

p o 2 'i 

v . = (8a,b,c) 

P r = maximum reflected overpressure (psi) 

P D — peak incident overpressure (psi) 

P a = ambient atmospheric pressure (psi) 
q Q = dynamic pressure (psi) 

P s = stagnation pressure (psi) 

The dynamic pressure referred to in Equations 8 a, b, and c is the result of the gross 
violent movement of air associated with the blast wind. It is smaller in magnitude than the 
blast wind itself when the incident overpressure is less than 70 psi, but is additive to the 
incident overpressure on vertical surfaces. Therefore, the stagnation pressure is equal to 
the peak incident overpressure plus the dynamic pressure, as shown above. 
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DDS 100-9 also lays out guidelines and empirical equations for calculating the beam 
resistance provided by a given plate-stiffener combination (stiffened panel). It uses the 
geometry of the plate and stiffeners to determine plastic section moduli at the mid-span and 
at the support. These moduli then combine with the dynamic yield stress of the material 
and the length of the stiffened panel to provide an overall resistance. The resistance of the 
overall stiffened panel is given as follows: 


^' a dy / \ 

r p = ~T~'( z ps +z pm) 

L (9) 

Rp = panel resistance (lbf) 

G dy = dynamic yield stress (psi) 

L = length of stiffener (in) 

Z ps = plastic section modulus at support (in 3 ) 

Zpm = plastic section modulus at mid-span (in 3 ) 

The yield stress used in this equation is a dynamic yield stress, including accounting 
for allowable plastic flow and dynamic load effects. It will be discussed farther in the next 
section of this report. This intrinsic resistance of the stiffened panel is then compared with 
the resistance necessary to withstand the desired blast loading. The required plastic load 
resistance for the blast is determined by applying the stagnation pressure (see Equation 
(8c)) over the breadth and length of the stiffened panel and dividing by an empirical force 
factor. This empirical factor is dependent on the natural period of the stiffener, the 


duration of the positive phase pressure pulse, and the ductility factor (p). The ductility 
factor is a dimensionless quantity determined in ship specifications that represents the 


amount of permanent set or plastic deformation that is allowable for various portions of 
the ship structure under blast loading. Equation 10 shows the formula for determining the 


normalizing force factor for the required plastic load resistance. 
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F 1 = empirical factor for normalizing plastic load resistance 
T n = natural period of stiffener (sec) 

T = positive phase duration (sec) 
p = ductility factor (dimensionless) 


( 10 ) 





The major design parameter for the stiffened panel is that its intrinsic resistance be greater 
than the required plastic load resistance for a blast designated in the ship design 
specifications. 

All of these design guidelines provide preliminary rules for current ship structural 
design. In addition to serving as a design tool to derive required plate thicknesses and 
stiffener moduli, the empirical relationships developed can also be used to evaluate existing 
structures. By reversing the design procedures, an existing panel can be analyzed to 
determine the blast pressure loading it can resist without failure. This method will be 
discussed further in following sections of this report. 

Beam Theory and Dynamic Loading 

Thus far, the methods of calculating the allowable loading on a stiffened panel have 
been empirical, derived from live fire data. To validate these results, an extensive study of 
the dynamic load limit of a stiffened panel was completed. The first step in this theoretical 
analysis was to define the maximum allowable stress in the panel before it is considered to 
have failed. As discussed, first yield is far too conservative a measure when blast loading is 
considered. First yield is defined as the elastic limit of a material, or the stress at which 
permanent deflection will result upon unloading. This point can be observed as the first 
knuckle in the stress-strain curve for steel shown in Figure 15. Due to the shape of a 
typical stress-strain curve for metals, a significant amount of plastic deformation can occur 
before the material approaches its fracture point. Strain hardening occurs during the first 
phase of this plastic deformation, such that the material is able to withstand still-increasing 
amounts of stress. The ultimate tensile strength of the material is defined as the highest 
stress the material can resist, corresponding to the highest point on the curve in Figure 15. 
Once this ultimate tensile strength is reached, the material continues to strain to fracture 
with no increase in applied stress. Since these results are obtained for an idealized 
specimen in a laboratory setting, using the ultimate tensile strength would be a bit risky as a 
measure of the fracture point of the material. Therefore, a quantity referred to as flow 
stress is used to indicate a representative stress for probable fracture. 






Conventional strain e 


Figure 15: Typical Stress-Strain Curve for Steel 

Flow stress equals the average between the yield stress and the ultimate tensile 
strength for ductile materials such as steel. It represents a stress high enough that 
continued increase in applied stress is likely to cause flowing plastic deformation and 
therefore lead to collapse or fracture of the material. The flow stress, a o , is used 
throughout the analysis presented here to represent the degree of stress that represents 
failure in a stiffened panel. 


In order to analyze panel geometry for plastic collapse, the bending moment 
representing the collapse condition must be defined. Jones (1997) defines the plastic 
collapse moment of any solid rectangular cross section as: 

M a-B-H 2 

4 ( 11 ) 
Mo = plastic collapse moment 
a = stress 

B = breadth of plate 
H = thickness of plate 

For the case examined here, the stress used in this formulation of plastic collapse moment 
is the flow stress. This moment is equated with the maximum internal moment across the 
breadth of the stiffened panel, a function of the loading. Solving for the load gives the 
maximum pressure the stiffened panel can support before plastic collapse. 
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Looking at a simplified geometry of the stiffened panel and treating it as a deep 
beam allow calculation of the internal moment as a function of loading. In two 
dimensions, the panel is modeled as a beam supported at the ends by large structural 
members such as girders. A single stiffener was added at the center of the beam to 
represent the supports on the panel. Actual panels encompass multiple stiffeners, but to 
simplify the mode shape and loading for this first estimate, a single stiffener case was used. 
A dishing mode shape was assumed for this simplified structure as shown in Figure 16 . 
Based on the assumed deflected shape of the structure, the center section of the beam was 
modeled to better define the boundary conditions. This modeled portion was cut at the 
zero-slope point of the deflected shape to allow modeling by a fully clamped end condition. 
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Figure 16: Simplified Stiffened Panel Mode Shape 

The modeled portion of the beam shown above was approximated as a two- 
dimensional beam with clamped end conditions to enforce the zero-slope deflection at this 
point. To properly represent the reaction forces on this modeled beam, the reaction force 
at the stiffener was assumed to carry the entire distributed load for the length of the 
modeled portion. The portions of the distributed load applied to the left and right of the 
modeled section of panel were assumed to be supported by the girders (or stiffeners) on 
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either end of the panel. This loading model represents a worst-case scenario where the 
stiffener is modeled as a large point load at the center of the beam. Figure 17 is a graphical 
representation of this model. 
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Figure 17: Modeled Portion of Stiffened Panel 

Superposition of the uniform load and point load cases from beam tables allow 
calculation of the total moment on the beam as a function of x-location. The development 
of this total moment is shown in Appendix E. The representative geometry and loading 
case used in the first section of Appendix E demonstrates the location of the maximum 
moment. As can be seen in Figure 18, the maximum moment occurs at the point of 
inflection produced (in both the deflection and the moment diagrams) at the point load 
representing the stiffener. For this simplified panel geometry, the maximum moment is 
located at the mid-span of the beam. Therefore, the exact moment at this point along the 
beam can be calculated and equated to the plastic collapse moment for the stiffened panel 
to solve for the maximum loading. 



Figure 18: Bending Moment in Modeled Stiffened Panel 
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The remaining two sections of Appendix E show the calculation of the maximum 
load supported by a stiffened panel as installed on the DDG-51 class and DD(X) class 
ships. Further discussion of the geometry of a representative stiffened panel for each 
platform is included in the next section of this report. Using the method described, the 
maximum uniform load, q, is determined for each geometry. Since the end goal of this 
calculation is a determination of the static pressure the panel can resist, this uniform load is 
divided by the length of the panel to re-introduce the third dimension to the problem and 
provide a result in terms of static pressure. 


This solution to the stiffened panel problem does not include any time-dependent 
or dynamic effects. Shock waves from blast loading propagate on a very short 
(microsecond) time scale. The positive pressure phase from the shock is only applied to 
the panel for approximately 2 x 10 4 seconds, and the peak pressure calculated in the section 
above represents the instantaneous maximum pressure reached during this positive phase. 
Therefore, the effects of dynamic loading must be considered when studying blast effects. 


The problem of the dynamic load condition can be modeled using some common 
strain-rate inclusion formulas such as the Cowper-Symonds relation (Equation 12). This 
relation outputs a “corrected” dynamic yield stress from input parameters including the 
static yield stress, the strain rate, and material constants. The strain to fracture or rupture of 
steel is used in this case to represent the desired maximum deflection for the panel. This 
strain is applied over the time period of the positive phase of shock pressure to obtain a 
strain rate. Empirical constants representing the material properties of steel are included in 
the relation. Strain rates typical for achieving the fracture strain in steel over the time 
period characteristic of a blast load will increase this dynamic yield stress by a factor of two. 



a odyn * a og* 


1 + 


(l *> 0 " 
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( 12 ) 


a odyn - dynamic yield stress 
G og = plastic flow stress 
£ dot = strain rate (= rupture strain/pulse time) 
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D - material constant (= 40 Hz for steel) 
n = material constant (= 5 for steel) 

The Cowper-Symonds relation is typically used to modify yield stress to include 
dynamic considerations. In this case, plastic flow stress is used in place of yield stress to 
represent the fracture condition imposed for this problem. It is unknown how these strain 
rates will affect the plastic flow stress of the material. It is possible that the steel will 
become effectively more brittle as a result of the high strain rate, therefore decreasing the 
fracture stress of the steel. Since calculations for plastic collapse moments are all based on 
the plastic flow stress, this is the important parameter to consider strain rate effects on. 
Plastic flow stress is calculated as a numerical average between the yield and fracture 
stresses for the material and therefore could likely remain close to constant under high 
strain-rate conditions. 

According to Brokenbrough and Johnson (1974), strain rate effects are also 
temperature dependent. For temperatures from —50°F to room temperature, increasing 
strain rate increases both yield strength and tensile strength in most steels, including high 
strength alloys and heat treated steels. However, at high temperatures (600°F) these effects 
are quite different. Increasing strain rate at high temperatures has very little effect on the 
initial yield strength of most steels. At the same high temperature, the ultimate tensile 
strength of all steels tested actually decreases with increasing strain rate. This is an 
interesting point to notice since many blast loads involve significant heat output as well as 
shock propagation. For the purposes of this study, however, it is assumed that the shock 
front propagates much faster than the heat from the blast can transfer to the steel. The 
thermal conductivity of air is approximately 0.024 Joules per meter-second-°Celcius and the 
thermal conductivity of steel is approximately 46 J/(m s °C). Based on these values, a 
minute amount of blast heat will have conducted from the air to the steel of the ship 
structure. Uniform heating effects such as those described above will not be a factor since 
heat would have to transfer to and through the steel material before these strain rate effects 
would be observed. Since the rate of heat transfer from air to steel is much slower than the 
shock front propagation, these high temperature effects can be neglected. Therefore, the 

assumption that both yield and tensile strengths increase with strain rate holds. 
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Plastic Collapse of Plates and Stiffened Panels 

To achieve more accurate results for the stiffened panel problem, Jones (1997) 
derives upper and lower bounds for a variety of panels and boundary conditions. The 
derivation of these bounds for a sample geometry is presented in Appendix F. The lower 
bound for a simply supported rectangular plate is calculated by deriving the bending 
moments in the x, y and cross directions and substituting these into Equation 13, the 
governing differential equation of equilibrium for the plate. Solving for p in this expression 
gives the static pressure to plastic collapse. 

-^M x + 2 • —— M xy + + p = 0 (13) 

dx dxdy dy 

The upper bound for the simply supported case is developed using a kinematically 
admissible transverse velocity profile. This profile, shown in Figure 19, assumes that 
regions I and II remain rigid, with plastic flow concentrated in the plastic hinges that form 
at the region boundaries. The admissibility of this profile is supported by observed failure 
modes in laboratory testing. 



Figure 19: Kinematically Admissible Transverse Velocity Profile 












An energy method is applied by deriving the internal energy dissipation in the plate, 
and the external work due to the uniform applied pressure. Equating these energies and 
solving obtain the upper bound for collapse pressure. This energy method is also used to 
solve for the collapse pressure in a partially restrained boundary condition case. To 
approximate the partially restrained boundary condition, a moment between zero and the 
plastic collapse moment, M 0 , is applied to the outer bound of the plate. This applied 
moment gives additional energy dissipation around the plate boundaries, which is added to 
the internal energy dissipation. This total energy is then equated to the external work 
energy (which remains the same from the simply supported case) and the plastic collapse 
pressure is obtained. 

Jones also looks at the dynamic plastic behavior of plates by including an inertia 
term in the governing equations of equilibrium. As a representative case, the static pressure 
derived above is applied as a uniform rectangular pulse of duration x. To simplify 
calculation, the plate is assumed to be square. The square geometry allows use of symmetry 
when assuming a kinematically admissible velocity profile. For cases where the applied 
peak pressure is greater than twice the static collapse pressure for the plate, the velocity 
profile must be divided into three distinct phases. The displacement solutions for each of 
these phases are added together to arrive at an overall displacement for the dynamic load 
case. The result of these derivations is presented in the first section of Appendix F. 

For the purposes of this study, this method of establishing net displacement based 
on a representative dynamic load is helpful to visualize the amount of permanent 
displacement a structure which does not rupture will undergo under these blast conditions. 
For example, using a plate thickness of 0.4375 inch, a one foot by one foot plate will deflect 
0.435 inch, or approximately one plate thickness. When the size of the simply supported 
plate is increased to more typical ship dimensions of 9 ft square, the resulting displacement 
is 7.62 feet. For a ship deck, this would represent a deflection spanning to within 2 feet of 
the overhead of the compartment. Some form of collapse or fracture failure of the 
structure would most certainly accompany deflections of this size. 




Further Considerations 

The theoretical methods presented above for solving the stiffened panel problem 
do not consider all factors influencing the final failure pressure of the panel. Each of "the 
methods used operates on a series of assumptions affecting the outcome achieved. These 
methods were chosen based on increasing complexity of the model used without resorting 
to numerical methods. It is recognized that many other factors must be considered to 
obtain truly accurate results. Often the affects of these conditions are only understood 
qualitatively and cannot be accurately modeled through empirical or simple first-principals 
based calculations. 

For example, the method of attachment of the stiffeners is likely to play a major 
role in this problem. Welding causes interesting modifications to the properties of the steel. 
In addition, most high-strength steels are worked to increase the yield strength of the steel, 
but in the process, the ductile zone of the member is reduced. This is especially true in the 
heat affected zone surrounding a weld. Therefore, although the stiffeners increase the 
overall strength of the panel until the onset of yield, the region of strain which the material 
can withstand after the onset of yield (usually about 8-1 Ox the yield strain for steel) can be 
significantly reduced due to the heat affected region. For example, steel ruptures when the 
strain reaches a fracture strain, s f , typically about 20% for bare steel. However, this strain is 
significantly less for material in a heat-affected zone. Paik et al (2003) show that actual 
weldment material has an even lower s f value, typically around 8%. 

DDG-51 and DD(X) Structural Analysis 

The question still remains; what is the level of mitigation necessary to make a blast 
mitigation system feasible? In order to determine a quantitative level of success for 
mitigation, real-world data for the strength of ship structures must be studied in 
comparison to a representative blast loading. Two candidate platforms were studied to 
determine their intrinsic blast resistance and the amount of mitigation necessary to prevent 
the spread of damage from a blast: the DDG-51 class destroyer and the DD(X) next- 
generation destroyer. 




The Arleigh Burke Class (DDG-51) destroyer represents the current state-of-the-art 
in Naval combatant design. Its structure includes a side shell and transverse bulkheads 
composed of high strength steel and a weather deck, longitudinal bulkheads, and stiffeners 
composed of ordinary or medium strength steel. The material properties for these steel 
components are listed in full in Appendix D. Analyzing the structure of this platform is 
very relevant since much of the impetus of designing blast mitigation systems results from 
the aftermath of the bombing of the USS Cole, an Arleigh Burke Class destroyer. 

DD(X) is the next generation of Naval destroyer, currendy in the detailed design 
phase. Although this ship does not yet exist in the operating fleet, it is a prime candidate 
for future installation of a blast mitigation system due to some of the ship’s unique design 
features. In the current design of DD(X), the side shell, weather deck, transverse bulkheads 
and internal frames are all composed of high strength steel with mild steel stiffeners. The 
longitudinal bulkheads, however, are composed of HY100 with HSLA-80 stiffeners. This 
stronger material is used as armor plating inboard of the Peripheral VLS cells to protect the 
ship from incoming blasts impacting the missile cells and causing collateral damage inboard 
of the launchers. Due to the advanced, tumblehome hull form used for this ship, the 
missile launchers are installed along the side hull rather than along centerline. The 
additional structure in the longitudinal bulkheads is both heavy and costly. However, there 
is space outboard of the missile launchers (due to the reverse curvature of the hull) which 
could potentially house a mitigation system. If blast pressures could be reduced outboard 
of the missile cells, a further layer of protection would be added to the hull. This 
protection could allow for a reduction in the weight and cost of the longitudinal bulkheads 
and reactive armor plating needed in the current design. 

For purposes of comparison, the structural arrangement of the side hull panels was 
analyzed for both the DDG and the DD(X). A variety of analyses akin to those presented 
above are calculated for each of these ships in Appendices E-I. Since most of the 
calculations presented do not specifically account for the support of the stiffeners in the 
panel, two geometries were considered: one using the panel as a whole and a second using a 





narrower panel spanning from frame to frame in width, but only from stiffener to stiffener 
in height. Figure 20 shows the relevant dimensions for each ship. 
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Figure 20: Stiffened Panel Structural Arrangement 
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Table 1: Stiffened Panel Properties 


Each of the methods of deriving collapse pressure for a stiffened panel described 
above was applied to the panels installed on these ships. In order to fully understand the 
geometry of these stiffened panels, the calculations were completed considering the panel 
to span from frame to frame horizontally and girder to girder vertically (LxD). From the 
dimensions shown in Table 1, it is observed that the geometry represents a nearly square 
panel for both ships. These calculations were then repeated considering the panel to span 









from frame to frame and from stiffener to stiffener (LxB). This panel has a much higher 
aspect ratio. The calculations involved in these analyses are presented in the corresponding 
appendices to this report. A summary of the findings of all calculations are presented in 
Table 2 and in Appendix I. 

When calculating the panel strength from girder to girder and frame to frame, the 
minimum pressure calculated is the result of the simply supported beam bending method. 
The maximum pressure from the first principles calculations is that resulting from the 
inclusion of a kinematically admissible velocity profile with a fully clamped boundary 
condition. The results from the empirical calculation methods used in the U.S. Navy design 
data sheets are two orders of magnitude greater than any of the first principle method 
results. Discovering the source of this discrepancy will be central to any future work on 
this problem. 

For the higher aspect ratio geometry between stiffeners, the beam bending method 
was calculated assuming the span to be twice the B dimension since the model for these 
calculations places one stiffener at mid-span and another at each end of the panel. Using 
the deep beam approximations where the span is significandy shorter than the depth, as is 
the case using the stiffener spacing, the bending moment calculations result in higher 
estimates of the resistance, especially when dynamic considerations are included. For this 
geometry, the values calculated for using the Design Data Sheets are still approximately an 
order of magnitude greater than those calculated using first principles. 
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Table 2: Comparison of Stiffened Panel Calculations 


For the purposes of this investigation, the peak overpressure applied to the panel in 
question is 130.3 psi or 8.985 x 10 3 Pa. The static equivalent pressure applied to the panel is 
312.4 psi or 2.154 x 10 6 Pa. This peak overpressure is below the failure pressure calculated 
using the U.S. Navy Design Data Sheets. The static equivalent pressure is less than the 
failure pressures calculated using the DDS’s for all cases other than the DDG-51 panel 
spanning from frame to frame and girder to girder. 
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The results of the stiffened panel problem provide a starting point to investigate the 
level of mitigation that can be achieved using water mist. The ideal situation would be to 
reduce the representative 500 lb TNT blast load on a ship bulkhead to a level below any of 
the pressures calculated above. This most conservative expectation is unlikely to be met by 
any mitigation system. However, achieving a level of mitigation which reduces the pressure 
loading to below the levels calculated using the dynamic considerations is a realistic goal for 
a feasible system. 





SCALING AND DESIGN 

Matching of Ship Failure Loads to Mitigation Level Achieved 

Based on the calculated load resistances of stiffened panels for two representative 
ships, a desired mitigation level can be identified. Evaluation of the feasibility of a water 
mist blast mitigation system for shipboard use begins with a comparison of the mitigation 
levels achieved through applying water mist to the mitigation necessary to prevent failure of 
ship structures. The same energy comparison method presented in Appendix A was also 
used to establish necessary mitigation levels for each of the shipboard panel geometries 
discussed. As a minimum, the failure pressure calculated using the fully clamped boundary 
condition and the energy method (Appendix F) were used to establish the amount of 
energy that must be absorbed from the impinging shock front. Results of these calculations 
are presented in Table 3 below. 


DDG-51 DDG-51 I3D(X) D1D(X) 

Units Full panel Between Full Panel Between 

Stiffeners Stiffeners 

A Pressure 

p /p % Rec 

A reduced/ A original ^ k O 

A Energy J 

Pa 2.087 xlO 6 1.479 x 10 6 2.037x10* 4.317x 10 3 

luctionof 96.9% 68.7% 94.6% 20% 

/erpressure 

/nr 1.435 x10 s 1.017 x 10 s 1401x10 s 2.966 x 10 4 


Table 3: Required Mitigation for Prevention of Ship Panel Failure 


From this table, it is observed that the mitigation necessary based on the full panel 
geometry is not possible with the levels of mitigation observed in laboratory experiments 
and numerical simulations. Even under the most liberal failure loading case calculated, the 
full panel geometry would require peak overpressures to be reduced to between 3 and 5% 
of their unmitigated values. This is about an order of magnitude greater mitigation than is 
achieved using water mist. It should be noted that these full panel calculations did not 
include the strengthening effects of the stiffeners in the panel and therefore are 


conservative. 





It is possible and quite feasible to mitigate the blast to the level, which prevents 
failure as calculated for the panel spanning from frame to frame and stiffener to stiffener. 
This stiffener to stiffener calculation is likely to be a closer approximation to the actual 
strength of the stiffened panel structure installed on U.S. Navy ships. Based on this 
consideration, the water mist blast mitigation system could be a useful defense to protect 
ships from blast loading on the side hull. However, the correct mist parameters must be 
designed to achieve the mitigation needed. In addition, the space used for the water mist 
system must provide adequate depth to callow the mist to exhibit its mitigation properties. 
A space providing only 1 or 2 meters of water mist prior to the shock front encountering 
the next bulkhead will actually enhance the peak overpressures of the shock front, creating 
a negative effect. This necessary depth, along with the proper mist properties, are essential 
to a feasible design. 

Mist Parameters 

It has been demonstrated that significant mitigation of the peak overpressures 
experienced in blast loading can be accomplished using water mist. Based on experimental 
results, it appears that peak overpressure can be reduced by 30% of its original value with a 
well-designed mitigation system using optimal droplet size and mass loading ratios. 
Numerical results for an optimized seeded gas indicate that pressure reductions as great as 
35% are possible. For a constant mass loading of 0.63 and optimized droplet size of 26-28 
pm, the maximum overpressure is reduced as shown in Figure 21. This figure shows an 
idealized pressure reduction, not accounting for the inherent decay of the shock pulse 
through air. The no particles case in this figure represents a pure vacuum, but does show 
the dramatic effect the water mist could theoretically have on peak overpressures due to 
shock. 
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Figure 21: Maximum Overpressure as a Function of Distance through Seeded 
Gas for M 0 = 1.49, r| = 0.63 (Schwer et al (2002)) 

One of the major design parameters in designing a mitigation system for shipboard 
applications is optimizing the mean droplet size of the water mist. For feasibility, the 
necessary nozzles to facilitate creation of the optimized droplet size should be readily 
available. Based on study of commercial water mist systems, the minimum practical mist 
size is about 15pm. The current U.S. Navy water mist system uses 100 pm droplets. From 
the experimental and numerical results presented in this report, the optimal droplet size 
falls between 20 and 40 pm. Therefore, the proposed mitigation system will use a mean 
droplet size of 25-30 pm. 


The mass loadings used in both the shock tube tests and the numerical simulations 
play a major role in the resulting mitigation. Some question exists as to whether a 
mitigation system should be designed for a specific mass loading or for a specific numerical 
particle density. For perfectly uniform droplets in a homogeneous mixture of mist, the 
difference between holding mass loading constant and holding number density of droplets 
constant would be identical. However, for a dynamic mist undergoing settling, introduction 
of new mist through the nozzles, the number of particles in a volume may not be constant 
for a constant mass loading and vice versa. 




Schwer et al (2002) investigated the effect of which of these parameters is held 
constant in numerical simulations. The results of these simulations are presented in Figure 
22 for a variety of droplet sizes. The plot on the left in this figure shows results for a 
constant mass loading, as many of the previous figures in this report have shown. The 
figure on the right shows the results from holding a constant number density of particles. 
The resulting mass loadings for each particle size are shown in the plot. These results play 
into the scaling of concentration used in this study based on the droplet size. The theory 
behind Schwer’s calculations was to neutralize the increased effect of vaporization for 
smaller droplets by reducing the total number of drops in the mixture. This method 
proved invalid since these vaporization effects are a function of surface area of the droplets 
to a greater extent than they are a function of total number of droplets. 
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Figure 22: Maximum Overpressure for Different Droplet Sizes with Mass Loading Held Constant 
(r| = 0.63) (left) and Number Density Held Constant (Schwer et al (2002)) 

Based on the results of the shock tube tests and numerical simulations presented, 
mitigation is most effective for mass loadings between 0.5 and 1.0. A large amount of data 
exists proving the mitigation properties for a representative Mach number for blast loads 
and a mass loading of 0.63. Therefore, this mass loading will be used for the proposed 
design, although it is recognized that this mass loading is predominantly a function of 
nozzle geometry and pressure on the system. For shipboard installations, optimization of 






the mass loading can only be achieved to an order of magnitude estimate. The actual mass 

/ • 

loading will be highly dependent on the current fire suppression system usage and pressure. 

Current U.S. Navy Water Mist Fire Suppression System 

The U.S. Navy has developed a machinery space water mist system that utilizes a 
modified high-pressure spray nozzle. According to the CPSMA (1997), The system is 
designed to produce high volumes of 100 mm droplet (mean diameter) sprays with very 
high spray momentum to achieve rapid suppression of large hydrocarbon pool or spray 
(Class B) fires. These nozzles emit 2 gpm at 1000 psi on a grid system. This system 
produces water flow application rates (0.06 — 0.07 gpm/ft 2 ) three to four times higher than 
commercially available systems. Although the mitigation effects are not observed to the 
same extent for droplet sizes in the 100 pm range used in the current U.S. Navy water mist 
system, modification to smaller nozzle apertures such as those used on many commercially 
available water mist systems could make such a system feasible for installation on ships. 

Input parameters into the design of a water mist blast mitigation system include the 
shape and size of the space to which the mist is applied and the pressure, flow application 
rate, and droplet size of the water mist itself. The 1000 psi system currently being tested for 
U.S. Navy platforms in a machinery room setting provides a good starting point for a 
feasible back fit design for future platforms. The mass loading (mass of particles/mass of 
air) for this system is on the order of 2.0, whereas the mass loading for the proposed 
mitigation system is 0.63. Therefore, the Navy system must be modified to reduce mass 
loading, a result that will be achieved simultaneously with reducing the mean droplet size 
applied to the space. 

Proposed DD(X) Design 

Based on the feasibility of the water mist fire suppression system already being used 
on Naval platforms to include a blast mitigation system, a design is proposed for the next 
generation destroyer, DD(X). DD(X) is a prime candidate for a blast mitigation system due 
to the geometry of the tumblehome hull form and the peripheral vertical launch system 
(PVLS). The sides of the hull slope inward from a maximum beam below the waterline. 
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Vertical launch missile cells are positioned directly inboard of the side hull in this design, 
creating a dead space outboard of the missiles due to the slope of the hull. Figure 23 shows 
the design of the PVLS modules, including the dead space outside of the launcher This 
dead space is approximately 2.5 meters deep in the athwartships direction at its deepest 
point, close to the waterline. 


Figure 23: PVLS Module and Outboard Hull Design 

The ship is designed with a bit more side-hull protection as a consequence of the 
location of the missile launchers. In order to protect the ship from damage due to 
unintentional detonation of the missiles housed in these launchers, current design uses 
reactive armor and heavy plating inboard of the launcher. This reactive armor is shown in 
Figure 24 along the back side of the module. This design is necessary to protect the ship, 
but cosdy in terms of space and weight. The proposed mitigation system provides an 
alternative means to partially reduce blast loading on the side hull of DD(X) could produce 
significant savings for the ship design. 




Figure 24: PVLS Schematic Diagram 








Figure 24 also shows the geometry of the space outboard of the missile launcher. 
This triangular space is the proposed location of the blast mitigation system. Potential 
concerns include the dimensions of this space being adequate to allow for mitigation of the 
impinging blast. The dimensions of the space at its deepest point fall close to the threshold 
where the water mist begins exhibiting mitigation effects. Therefore, a close look at the 
exact dimension necessary on a large scale application to achieve mitigation effects is 
essential to the continuation of this design. Unfortunately, enhancing this issue, the top 
portion of this space does not provide the same depth. In terms of a system to protect 
against small boat threats, this is not a major issue. The top of this space is well above the 
waterline of the ship and outside the reach of an impinging blast wave, but depth of the 
proposed installation is still a major concern for future work. 

This blast mitigation system is applicable for back fit or installation on any Naval 
platform using the water mist fire suppression system. DD(X) is proposed here due to the 
inherent dead space just inside the hull, but since the water mist blast mitigation system 
only requires the installation of additional nozzles in the spaces direcdy inboard of the hull, 
it could be installed on any platform. 

The premise of operation for such a blast mitigation system would be to light off 
the mist nozzles in the spaces inboard of the hull when an incoming threat was eminent. 
This could mean an incoming missile threat in an at-sea scenario or an approaching small 
boat in a littoral or in port scenario. Constant operation of the nozzles would not be 
necessary; further investigation would be required to define a cycling process to maintain 
the necessary mass loading of mist in the air. 







CONCLUSIONS AND RECOMMENDATIONS 


The recent trend of using fine water mist systems to replace the legacy HALON- 
1301 fire suppression systems warrants further study into other applications of the water 
mist systems. Shock tube testing and numerical simulations indicate that fine mists (20-30 
pm droplet size) may reduce peak overpressures of a shock wave traveling through a space 
by up to 70%. Such pressure reductions could be used to mitigate the destructive effects of 

a shock wave (initiated by an explosive device) traveling through a structure. Although the 

/ 

mitigation effects are not observed to the same extent for droplet sizes in the 100 |tm range 
used in the current U.S. Navy water mist system, modification to smaller nozzle apertures 
such as those used on many commercially available water mist systems could make such a 
system feasible for installation on ships. 

Currently these blast mitigation effects have only been demonstrated in small-scale 
shock tube tests and computer simulations. Uncertainty exists as to the scalability of such a 
system. The intention of this research was to fully investigate the applicability of such a 
blast mitigation system for shipboard use. A study of the degree of mitigation necessary to 
make a system practical for shipboard installation was conducted. In addition, a theoretical 
study of the mechanisms of blast mitigation using water mists was undertaken to explain 
the phenomena and possible methods of scaling and application. Current research into 
water mist blast mitigation was thoroughly reviewed and examined for trends to aid in 
guiding the direction of future work. Based on these studies, it was determined that water 
mist systems do provide enough peak blast pressure reduction to warrant possible 
installation on warships. 

These trends provide the necessary information to design a scalable blast mitigation 
system with applicability to real-world shipboard spaces. Given the recent trend toward 
tumblehome hull forms in future Naval Combatant designs, there exists strong applicability 
of this system in the “dead” spaces created by the shaping of the tumblehome hull. The 
end goal was a feasible design of a blast mitigation system to be used in the outermost 










spaces of Naval Combatants to protect interior vital system spaces. However, concerns 
over the depth of mist that a shock wave must travel through in order to achieve mitigation 
limit the applicability of this system based on current data. 

A blast mitigation system for installation on U.S. Navy ships using water mist is a 
still a viable future technology. Based on the strength characteristics of a current and a 
future destroyer, the mitigation system would be required to reduce peak overpressures 
from a blast by between 30% and 80% of the initial blast loading pressure. Depending on 
the strength of the stiffened panels in the ship, this mitigation could reduce the loading on 
the ship structure to a level that makes the probability of plastic collapse or failure of the 
panel minimal. More research into the applicability of scaling the current laboratory and 
numerical results to a shipboard space are needed to continue design of the system. 

Recommendations for Future Work 

Further development of a water mist mitigation system would require quantitative 
mitigation results for the proposed design using shock tube testing on an actual water mist. 
All laboratory data presented here was for some form of dry particle-air mixture. The 
results for water droplets are all based on numerical simulations validated against the dry 
particle results. A method of suspending a representative water mist in the driven section 
of a shock tube must be developed to accomplish this testing. Outa et al (1976) discuss the 
difficulties involved in developing a technique for injecting particles to achieve a 
satisfactorily uniform dispersion in the gas during shock tube testing. These conditions are 
amplified by the requirement to create the water mist and ensure that the mean droplet 
diameter of the mist meets the design specifications for the mitigation system. Design, and 
testing of a method to run shock tube simulations on water mists of varying concentrations 
and droplet diameters should be accomplished. The results of such testing would validate 
numerical simulations of the mitigation effects and allow more detailed design of a full-scale 
mitigation system. 

Scaling of these shock tube tests to a more typical shipboard space must be 
investigated to further validate current numerical models. This could be accomplished by 






running live fire tests with water mist in a space and recording pressure profiles along 
various locations inside the space. This data would be invaluable to determining the spatial 
effects of using the water mist in a space that has a breadth and height comparable to its 
depth. 


In addition, a close study must be conducted to quantify the exact depth of space 
required to achieve mitigation effects from a water mist using the parameters optimized for 
the proposed system in this report. Actual shipboard spaces should be used to model the 
scaling effects of moving from a small shock tube to a large (in terms of breadth and 
height) space. The important parameter to achieve out of these tests is the necessary depth 
required to ensure that the mist provides adequate mitigation effects and does not cause 
enhancement of peak overpressures under any conditions. Determining this depth 
parameter will define whether the mitigation system is feasible for future installation. 

Future work in the problem of blast loading and mitigation systems for ships 
should also include more analysis into the stiffened panel problem. Based on the 
theoretical calculations presented in this report, the loading to failure of a typical ship 
stiffened panel is on the order of 10° Pascals. However, the failure loading calculated using 
the empirical formulations in U.S. Navy Design Data Sheets is about two orders of 
magnitude greater than this. Although the theoretical calculations used in this analysis 
involve many simplifying assumptions including mode shape and velocity profile, a wide 
variety of methods were studied. Of the parameters varied among these methods, 
boundary conditions had the largest effect on the resulting pressure load to failure. None 
of these parameters or assumptions account for a discrepancy as large as that between the 
theoretical results and the empirical results. The cause of this discrepancy should be 
identified and understood before lending further credence to the results for the actual 
strength of stiffened panels under dynamic blast loading. 

Finally, for further development of the specific design proposed for DD(X), the 
stiffened panel analysis presented should be repeated using the scantlings of the armor 
plated bulkhead located immediately behind the PVLS. The specific geometry of the space 





outboard of the PVLS launcher should be entered into numerical simulation tools to 
establish a mitigation pattern for this proposed installation. Based on the mitigation levels 
likely using a water mist system, the longitudinal bulkhead could be optimized. By 
optimizing this bulkhead, an analysis could be conducted on the difference in ship 
displacement, cost, and fabricability based on the installation of the blast mitigation system. 
A significant savings in terms of ship design impact would provide the true determination 
of the feasibility of the water mist blast mitigation system. 
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APPENDIX A: PROPERTIES OF WATER MIST 


Properties of Water Mist: 

ps/p is normally - 10 A 3. For fresh water mist in room temperature air: 


P water - 3 Pair'” 1*293- ^ 

m m 


P water 


= 773.395 


^water 


HeatofVap := 539-^- 
gm 


C air :=0.25- 


^ratio 4 


HeatofVap = 2.257x 10 — 
kg 


Mass loading for current U.S. Navy water mist fire suppression is as follows: 


FlowAppRate := 0.065- 


* A 2 

min-ft 


MassAppRate := FlowAppRate - p 


FlowAppRate = 2.648- 


. 2 
min-m 


MassAppRate = 0.044- 

2 

m s 


Remove the time dependency by assuming a settling rate of 1 meter per minute. 


Concentration := MassAppRate -1.0-- 


To vaporize this amount of mist: 


Concentration = 2.648— 
3 
m 


Energy := HeatofVap-Concentration Energy = 5.977x 10 — 


This assumes the water mist is already at 100C. To get from room temperature, must calculate heat 
transfer to droplets: 


k a := 0.024- 


Nu(D) :=h aiV — 


Kvater 0.586 


hair * * ^ 


hr-ft 2 -R 


h air = 102.209 


a(25 10 6 -m) 


Nu„:=-0.10: Nui nr> := 0.421 
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For a representative temperature shift, use AT = 80C 


Qcdot( D ) ;= Nu(D)-7t-D ^ vater A T 


t (25 10 6 -m) 


AT := 80-K 


= 3.943 x 10 4 W 


per drop 


Now need to know how many drops are in 1 m A 3: 


M ass dr 0 p(D). Pwater' - 


-(25 10 6 *m) = 8.181 x 10 12 kg 


n drops ^ 


Concentration 

Mass drop( D ) 


.(25-10 6 -m) 


= 3.237 x 10 


Power temp (D) : = "drops( D ) Qcdot( D ) 


Power temp (25-10 6 -m) = 1.277x 10 8 — 


The time period for transfer of this heat to the droplets remains to be known. This time 
dependence is the primary reason that it is likely that vaporization is not a major contributing factor 
to the overall energy absorption. However, for example, assume the heat transfer does occur over 
the time of the shock pulse: 


t := 0.2-10 3 -s 


Energy t „ mn (D) := x-Power temri (D) 


Energy tem p(25-10 ^-m) = 2.553 x 10 4 — 


Energytotal (°) := Energyvap + Ener gytemp( D ) 


Ener gy total 


(25-10 6 -m) 


= 6.002 x 10 6 — 


For comparison, calculate the specific energy in the representative pressure pulse impinging on this 
volume: , x 

-(‘-‘do) 


Po := 130.3-psi 0o := 772.510 sec t Hn := 0.029sec 


Pi(t) := Po e 


Area := 1-m 


Volume := 1-m 


Speed_of_sound := 344-— 
s 


Energy v 


Pj(t) ■ Speed_of_sound dt 


0.029sec 


Energy =2.387x 10 


Energy m j. t := 6.002-10 - 


NetEnergy := Energy wave -Area - Energy mist -Volumt 


NetEnergy = -5.763 x 10 J 
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Shock energy in 




# 


m 


m 


This method shows that the energy of vaporization of the water mist is over an order of magnitude 
greater than the energy of the pressure wave. For verification, use the static equivalent pressure from 
DDS 100-9 in a rectangular pulse form to recalculate wave energy: 


P a := 14.7-psi Pj := Po 


Pf := 2'Pp 


' 7-P a + 4-P^ 
v 7'^a + ^1 ) 


P r = 697.429psi 
P r = 4.809 x 10 6 Pa 


% : = 

% = 


V 7-P a+ P 


182.012psi 


.255 x 10 Pa 


i) 


P s : ( P 1 + %) 

P s =312.312psi 
P s = 2.153x 10 6 Pa 


5 J 

Energy2 wave := P s *t ^peed_of_sound Energy2 wave = 1.481 x 10 — 

m 

This wave energy is slightly lower than that calculated above for the exponential decay, but of the same 
order, validating the results above. 
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Recognizing that the mass loading for the current Navy system is much greater than most 
commercial systems using smaller droplet sizes, recalculate based on a smaller concentration of 
water mist which is more representative of the mist applied through the smaller aperture nozzles 
needed for this droplet size (reduce concentration to 1/4 since reducing mist diameter by a factor 
of 4): 

„ Concentration 

Concen :=- 

4 


To vaporize this amount of mist: 


Energy ya p 2 := HeatofVapConcen 


Energy vap 2 = L494 * 10 6 

m 


This assumes the water mist is already at 100C. Based on results of previous calculations, 
neglect the effect of raising the temperature of the drops from room temperature. Even at this 
reduced mass loading, the specific energy of vaporization of the mist is an order of magnitude 
greater than the energy of the impinging wave. 





APPENDIX B: DYNAMICS OF GAS-PARTICLE 

MIXTURES 


From Marble (1970): 

Ideal Gas Properties: Particle Cloud: 

velocity uj particle radius 0 

temperature t velocity u p j 

density p temperature T pi 

particle mass n: 

number particles/unit volume n m 

density p p = n p • m 

Due to non-interaction, all particles in a local volume have the same velocity vector and 
temperature - individual variables could deviate randomly because of Brownian motion 
(particles too large) or random initial conditions (damped early in the history of the problem) or 
due to multiple particles sizes in mix (assume uniform here). 


Continuity: 

(volume of particles considered negligible) 


0 0/ \ r 

-p + —(p ■ uij = C 

0t 0X[ 

“Pp +— (pp * Upi) = 0 
0t 0?q 


By motion equilibruim for local gas motion/acceleration, during the coupled motion of the 
gas and particle cloud, the particles adjust to the slip values (difference between gas velocity 
and particle velocity) corresponding to the local gas acceleration, while traveling a distance 
small compared to the characteristic dimensions of the problem. (Similar formulation for 
thermal equilibrium). 


Force exerted by a single particle moving through the gas: 
Force per unit volume with n non-interacting particles: 


where 


6 • 71 • 0 • \i 


fp = 6 71 • CT • n • (u p - u) 

Fp = n • 6 • 7i • a • |x ■ (up - u) 
( U P " u ) 


^P “ P P ' 


is the characteristic time, or the time required by a particle to reduce its 
velocity relative to the gas (slip velocity) by e^-1 of its original value 
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Limiting form of system results as t goes to 0. In this case, slip velocity (up-u) also goes to 0, but 
the volumetric particle force remains finite. This limit results in: 


U pi = Uj 


D (~] 

and therefore --— = Constant 

Dt 


In words, this means that particles are "fixed” to their initial mass of gas. If the initial particle 
distribution is uniform, pp/p = K = constant. Equations of Motion become: 


0 0 —10 1 
- Ui + u j * — u i = 7"-x-P + E-7P 

at a^ U + k )p 0xj U + K ) 







Under this model, the system behaves like a perfect gas with modified properties: 
Density becomes: p' = (l + k)p 


Kinematic Viscosity: 
Reynolds Number: 


, _ H _ v 

0 + k) P (l + k) 

u 0 ■ L (l + k) • u 0 • L 

Re' =-=- 

v* v 


Mach Number: 



U 0 V 

(l + k)| 

l c v ) 

a 


K • C 


1 + — 

> 


c p 




Effective flow operates under larger Reynolds and Mach numbers 

Shock structure consists of complete equilibrium of gas and particles ahead of and behind 
the normal shock. A zone of gas-dynamic shock and a zone in which velocity and 
temperature equilibration occur compose the shock structure itself. To separate these 
zones, the velocity and temperature relaxation times must be long compared with the 
molecular collision time of the gas. 
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e 


^.(Tj >1 
U 


ps = density of the substance constituting the particles 
p = density of gas 
a = particle radius 
1 = mean free path 

ps/p is normally ~ 10 A 3. For fresh water mist in room temperature air: 

Jr p 

Ps := 1000- — 


m 


p := 1.293 — 
3 
m 


Ps 

— = 773.395 


Marble postulates that the gas-dynamic shock structure is unaffected by the particle cloud 
even though the shock may be thin relative to the particle radius. Because the volume 
fraction of the particles is negligible, at any time less than 1% of the shock front is being 
punctured by particles. 

Sommerfield extends Marble’s postulations using numerical calculations and experimental 
results. He first postulates based on loading ratios of the gas-particle mixture for both mass 
and specific heat: 


Qvater 1-00* 


cal 

gm ■ K 


C a ir 0.25 - 


cal 


gm • K 


Q-atio :_ 


Qvater 


Qatio 


= 4 


Mass loading for current U.S. Navy water mist fire suppression is as follows: 
FlowAppRate := 0.065* ——- FlowAppRate = 4.414x 10 ^ms ^ 

■ a 2 

min ■ ft 

-2 -1 

MassAppRate := FlowAppRate • p s MassAppRate = 0.044kgm s 

Remove the time dependency by assuming a settling rate of 1 meter per minute. 


Concentration := MassAppRate *1.0* - Concentration = 2.648kg m^ 

m 

Therefore, the mass loading (mp/mg)is: 

_ _ _ Concentration 

MassLoading :=- ; MassLoading = 2.048 


P 

Therefore, the current Navy system has a mass loading higher than that for which the studies at hand 
have calculated mitigation effects. 
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Recognizing that the mass loading for the current Navy system is much greater than most 
commercial systems using smaller droplet sizes, recalculate based on a smaller concentration 
of water mist which is more representative of the mist applied through the smaller apeture 
nozzles needed for this droplet size (reduce concentration to 1/4 since reducing mist 
diameter by a factor of 4): 


_ Concentration 

Concen :=- 


MassLoad 


Concen 


MassLoad =0.512 


4 


P 
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td is the time delay for the pressure wave arrival at distance R 
Ro is the shortest radial distance 
c is the shock wave velocity 


R- R 0 

l d( R ’ R 0’ c ) := “ 

t d ( 10,0,344) =0.029 

Can plug td into P(t) to get actual shock front pressure - (t-t di ) 

P oi := 4.82910 7 0i:= 772.510" 6 t di := 0.02' P(t):=P oi -e 01 



From General Specifications and DDS 100-7f, Missile Blast Pressure Equivalent: 


P = 



0.0225 ^ 
sin(ct) ) 


A 


P = static equivalent pressure (psi) 

T = thrust of the missile (lbs) 
a = angle of incidence (degrees) 

A = impingement are of the surface (in A 2) 






From DDS 100-9, Nuclear airblast design for surface ships: 


P a := 14.7-psi Pj = 130.316psi 

f 7-P a +4-PO A Pi 2 ^ 

r ‘ = ' V { 7-Pa+Pl ) q °''2\7.P a + P l J 


P s ( P 1 + ^o) 


P r = 697.54 lpsi q 0 = 182.045psi P s = 312.362psi 

P r = 4.809x 10 6 Pa q 0 = 1.255x 10 6 Pa P s = 2.154x 10 6 Pa 

From Remennikov, A. (2004): 

Empirical Method: 

Radial Stand-off Distance (m): r := ic Charge Weight (kg): w o := 50 C 



These static and reflected overpressures are close to those calculated by Nonaka 
above, validating that method of determining peak overpressure. The results using 
Ramajeyathilagam's empirical formula are over two orders of magnitude higher, although the 
general form of the equation is correct. An error is suspected in the empirical constant used in 
this formulation. 

As expected, the reflected overpressure from the graph is significantly lower than 
that calculated from DDS 100-9. This is due to the open, single face geometry used in the 
empirical method, vice the confined space. Based on these calculations, a representative blast, 
roughly equivalent to that experienced by the USS COLE (~500 lbs TNT) is as a comparison for 
the level of loading ship structures can resist without failure. 
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Reflected Overpressure 


Dynamic Pressure 


Static Equivalent Overpressure 






APPENDIX D: MATERIAL CONSTANTS 


Properties of Steels Used in Ship Structures 


For Medium Steel (MS or OS) Components: 


a yms := 34-ksi o yms = 2.344x 10 Pa 


°ten ms := 58 ksi 


E rupms ' 


For High Strength Steel (HSS) Components: 


a yhs := 51 ‘ ksi 


°ten hs : = 7 °- ksi 


o yh$ = 3.516x 10 Pa 


E ruphs 


:= 0.23 


For HY*80 Components: 


a y80 := 80-ksi a y80 = 5*516x 10 P a 


For HY-lOO Components: 
J yl00 := °y!00 : 


G ten 100 


:= 110-ksi 


E rupl00 


:= 0.18 


v := .30 


E ms := 29.610 ksi 


Pms : = TO ^ 

m 


E hs :=29.610 3 ksi 

Phs :=7833^ 
m 


E g0 := 29.610 -ksi 


J ten_80 := 95 ksi e rap80 := °- 20 P 80 := 7833 


) 3 ! • 

kg 


E 100 := 29.610 -ksi 


P 100 := 783 ^ 
m 


E^ = 2.041 x 10 U Pa 

r _ Ems 
ms ' 2-(l + v) 

E^ s = 2.041 x 10 H Pa 


Ghs : = 


G hs 


2 (l + v) 


11 . 


Eg 0 = 2.041x10 Pa 


G 80 :_ 


-80 


2-(l + v) 


E 100 = 2.041 x 10 H Pa 


G 100 


-100 


2-(l + v) 
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Panel Characteristics for DDG-51 and DDOQ 


Panel Geometry 




Deck or Girder 


p 

Frame 

Stiffeners 

Frame 


Deck or Girder 

1 

___ 



1 


L 

— 


Deck 

Height 


DDG-51 

T g := .3438 in 
B g := 27.39in 
L g := 8-ft 
H g := 9.13-ft 
NoStiff g 3 
A g :== 1.81- in 2 
BFg := 3.94in 
I g := 6.5in 4 


Panel Characteristics 
T = Plate thickness 
B = Stiffener spacing 
L = Frame spacing 
H = Deck height 
Number of Stiffeners 
Area of Stiffener 
Breadth of Flange 


DD(X) 

T x ; = .4375 in 
B x := 23.99 in 
:= 7.546ft 
H x := 9.99 ft 
NoStif^ := 4 
A x := 2.11-in 2 
BF X := 3.94 in 
I x := 14.2-in 4 


Moment of Inertia 









DDG-51 


Panel Stiffener Characteristics 


DD(X) 


D g := 5.90-in 
TW g := .17ain 


TFg := .215 in 


A w t bg»TW g (D e -TF g ) 


A flag : TF g“ BF g 


Depth of Web 
Thickness of Web 
Thickness of Flange 
Area of Web 
Area of Flange 


D x := 7.89 in 


TW X := .170 in 


TF X := .205 in 


A wcbx : = TW x -(D x -T F x ) 


^flax : TF x BF > 




( a yhs + a ten_hs ) 
2 


a og = 6-05 x 10 psi 


Plastic Flow Stress 


( a yhs + a ten_hs ) 
° ox ' 2 
4 

a ox = 6 05 x 10 psi 


a 0 g = 6.05 x 10 psi 


cr ox = 6.05x 10 psi 
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APPENDIX E: PLATE BENDING CALCULATIONS 

Calculation of Bending Moment for a Stiffened Panel using Sample Geometn 


Sample Plate Characteristics: 
Thickness: j^ : = .4375 in 


k'-= 12-in 


B := 12* in 


7 yhs + a ten_hs ) Plastic Flow Stress 


Sample Blast Characteristics: 


p := 130.316psi 


p = 8.985x 10 Pa 


t:= 0.2-10 *s 


Plastic Collapse Moment: 

Plastic collapse moment of solid cross-section: 


Geometry: 

Assuming a dishing mode shape between stiffeners and frames: 

Undeflected Shape 



M 0 = 3.925x 10' N-m 


Loading 



Deflected Shape 



Modeled Portion 







1 



Approximate the modeled portion of the stiffened panel as a beam (2D) with a fixed end condition to 
enforce the zero-slope inflection point mid-way between stiffener and frame: 


Mi 


left 





Li_Jz_Jz_iLi_Jz_Jz_Jz_Jz_ J 


n 

rrr^ 

- > X 



Stiffener P\q*B 


right 


Assume that reaction force for the stiffener represents the entire distributed load, q (= pressure per 
unit width) * B (length between stiffener and next structural member). 

q:=p-L q=2.739xl0 5 — L = 0.305m 

2 

s 

From Beam tables, superposition of two loading cases to get the total moment: 


M uniform( x ) : =^-( 6B ^-B 2 -6x 2 ) 
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^total( x ) : ^uniform^ x ) + ^concentrated ( x ) 


^center * ^total^J 
^ends ^total ( B ) 



M center = ~212x 10 3 N-n, 
Mends = 106x 


This representative loading case is used to show that the maximum moment will be at the at the 
point of inflection produced at the stiffener (x = B/2). Now, symbolically solving for this moment 
and setting equal to the plastic collapse moment, Mo: 

M 0 = 3.925 x 10 3 N.m 

2 2 

^unif ~ 24 M conc “ g M tot ” M unif + M conc 


M tot = M Q For Plastic Flow / Failure 


M 


tot 



Neglect negative sign - Mo is for either pos or neg moment 






Calculation for Panel from Frame to Frame and Deck Gitdet to Deck Girder 


Plastic Collapse Moment: 


a og B g'^g 4 

M og := 4 M og = 1.939x 10 N-m 


w °ox^x^x 4 

M w :=- IVL V = 2.962x 10 N-m 


ox ■ 


M, 


unif = 


q 0 B 

24 


M, 


-q 0 - B 


cone 


^tot ^unif + ^conc 


M tot = M c For Plastic Flow / Failure 


-12 

M tot = —-q 0 -B Neglect negative sign - Mo is for either pos or neg moment 


q og : 


M og ‘12 


B„ 


4 og 


p og - Y 


q og = 3.005x 10 4 — 
s 


P 0g = 1.232X 10 Pa 


q ox * 


Mqx'12 


Pox • 


B x 

^ox 


q ox = 3.833 x 10 4 — 


p ox = 1.667x 10 Pa 


Now, determine the need for inclusion of dynamic terms by looking at kinetic energy 


Average strain rate to rupture will be: 

£ ruphs 


t:= 0.2-10 3 -s 


£ dot • 


£ dot = h!5x 10 Hz 


For strain rate effects, include Cowper-Symonds relation: 

For steel: ^ , 1 8 

D := 40 — n := 5 °yhs = 3.516x 10 Pa 


a ten hs = 4-826 x 10 Pa 


c odyn : °og 


1 + 


r«.dotV 

{—) j 


a og = 4.171 x 10 Pa 


^odyn = 1 234 x 10 Pa 


°odyn = l-789x 10 5 psi 
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a ogd 


^odyn'^g'^g 


a ogd 


A oxd 


a odyn' L x‘ H x 


%gd - = 


p ogd 


^ogd 

B g 2 

%gd 

V 


q ogd = 8.887x 10*^ 


Pogd 


= 3.645x 10 Pa 


q oxd *” ‘ 


M oxd- 12 


%xd 


%xd 

Poxd *” Poxd 


oxd =8.76x 10 4 N- 

= 1.134x 10 5 — 

2 

s 

= 4.929x 10 4 Pa 
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„, a odyn ^g'^g 

:=- M ogd = 5.735x 10 N-m M nyH :=--— : - = 8.76x 10’N-m 


A ogd ■ 


A^r ., a odyn ,L x'^x 

4 oxd := -:- 


A oxd 


%gd • 


Pogd : 


M ogd 12 


B_ 


^ogd 

"v 


^ogd 


Pogd 


5 kg 


= 3.555x Iff 


= 1.458x 10 Pa 


%xd * 


M oxd* 12 


Poxd * 


B x 

^oxd 


q ox d= 7.078x 10 5 ^ 


Poxd = 3.077x 10 Pa 







APPENDIX F: PLATE PLASTIC COLLAPSE 
CALCULATION 


Dynamically Loaded Plate Analysis from Jones (1997) for Sample Geometry 


Plate Characteristics: 

Thickness: 0.4375 • in 




B := 12 in 


A ,2 ,2 ,4 

^ 4 w(x 5 y)+ 2 ■ ~ 2~2 w (x>y) + ^4 w(x,y)= -J 
dx dx dy dy D 


D := ■ 


E-H 




12 


G :=- - -- 

2 • (l + v) 




( c y + c ten) 


Plastic Flow Stress 


Blast Characteristics: p o :=130psi t:=0.210 3 s 

Static Plastic Collapse Pressure of Rectangular Plates 


Plastic collapse moment of solid cross-section: 


M 0 := 


Jo-BH 2 


Plastic collapse moment of solid x-section of unit width: 


Simple Supports: 




JoH 2 


Use flow stress, cro, rather 
than yield stress to get 
plastic bending moment 


M 0 = 1.288x 10 4 kgms ' 2 



; 

X 

- -- -.y 

\ 

f 

y 

s£---i 


2B 


2L 


rn 




^ ^ 


j/ X X J/ \i/ si J/ X X \l J/ I 


XX 


T 


XX 


H 
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Angular velocity at the outer boundaries (simply supported) are: 

±9 Wdot 

dt B • tan(<|)) 


Internal Energy Dissipation in plate: 


^dot ~ 4 * M 0 • W ^ g + c °t( < l ) ) j 


External work due to uniform pressure p_upper: Ej 0 t = 2 • B Z * Wj ot • Pupper 


L tan(<|)) ^ 
B 3 ) 


3 : = — (aspect ratio of plate) 

L 

Equating these energies and solving for pjjpper: p U pper = 6 • M 0 ■ 


( 1+ _eJ> 

B 2 • (3 - p • tan( 4 >)) 


Assumptions are made to solve for the angle of the plastic hinge line along which the velocity profile 
follows. Namely, the minimum p_upper occurs when the derivative of p_upper with respect to tancj) is 
zero. Solving for tan<|), get: 


i(*) = -P + V3+P 


which can then be substituted into the energy solution for pjjpper 


Pupperss ‘ 


Pupperss — 8.317x 10 Pa 


B 2 • vv 3 + ( 3 2 - P, 


Partially restrained supports (applying moment m to outer bound): 


0 < m < M 0 


m_supp determines the level of restraint on supports: 
m=0 for simply supported, m=Mo for fully clamped 


n^supp •“ M 0 


Gives additional internal energy dissapation around plate boundaries: 

n_, - f » . k\ 


Bdot b = 4 • m • Wj 0 t • 


B*tan((|)) B J 


Add this to internal energy for the simple support case (from above) for: 
Ddot tot = 4- M 0 - W dot • f-t + -L.) • f 1 + ^ 

V B tanH); ^ j 
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The external work from p remains the same as above. Equating and solving for p: 


Pupper_fc 6 ' M t 


H^SUpp ^ 

~m7 ) 


B 2 .y 3 +p 2 -p 


6 - 1-2 

Pupper fc = 1-663x 10 kgm s 


For the fully clamped case, the lower bound is: 


Plower fc4 • M 0 


(, + P 2 ) 


Plower fc~ 1.109x 10 kgm s 


For a fully clamped square plate, p goes to 1 and a numerical solution for the exact static collapse press 
can be found: >, 

Mo 6-1-2 

Pc fc*“ 10.71- Pc fc= 1-485 x 10 kgm s 

,L 2 

Dynamic Plastic Behavior: Governing equations for rectangular plates 

p, r\ ri 

—Q x + — Q y +p = p • —^ w note inclusion of inertia term in transverse direction 

Sx 3y di 


— M x + —M xy - Q x = 0 
0x dy 

— My +-M xy -Qy = 0 
0y 9x 


Curvature Definitions in x and y directions: 


K x = -2 w K y =-2 W 


First look at the case of a rectangular pulse of pressure p_o lasting from time=0 until 
time=x 


t := 2.0 • 10 s 


p(t) := p 0 - dunif(t,0,x) 



M0 4 2-10 4 





Simply supported square plate: Static Collapse Pressure 

6M 0 5 _i 

Pc_ss : = —— Pc_ss = 8 - 317 x 1 0 kgm 

T 




(x+ y) • m 

V2-L 


(from above) 


— w = W dot • (1 - z) 

at 


Plastic hinges form symmetrically along diagonal > 
and y axes, therefore can consider a single ' 
quadrant. 







These equations are plugged in for each time period (0-t, t-T) where x is the end of the pressure pulse 
and T is the end of motion in the plate. 

Po 

h SS • Tfss — 1-078 

Pcss 

For values of r\ greater than 2, a yield violation develops near the plate center, requiring a modification of the 
above transverse velocity field into three distinct phases: 



w dot = W dot 0<z ^ 


w dot - W do t • 


(1-z) 

0-Co) 


4o<z< 1 




(solution admissable for r| greater than 2): 

■ Tl SS _ 2) 

4-p 


Wf2 ss hss ' Pc ss ' T 


( 3 


w f2 ss = 0.011m 
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Lower bound (Johansen Yield Condition): 


( 


Mx(x):=M 0 


V 


hr') 


( 


My(y) :=M 0 


1 - 


V 


2 ^ 

y 1 

Bpr 2 J 


Define: 




M xy (x,y) := 


-M 0 - x-y 
Bpr'Lpr 


Governing Equation: 


-^M x + 2--—M xy + -^M y + p = C 
Sx 9x6y dy 


Substituting and solving for p: 


Plower ss 2 ' Mo - 


, L f L V 
1 + — + — 

B UJ 


4 -1 -2 

Plower ss ~ 2,829x 10 kgm s 


Upper bound (using energy method and kinematically admissable transverse velocity profile): 


Velocity Profile: 


Region I (triangles at ends) w dot = W dot 


Region II (center section) 


w dot = W dot - 


(Bp r • tan 4) - Xp r ) 
Bpr • tan((j>) 

(Bpr ~ y) 

Bpr 


Each of these regions remains rigid with plastic flow concentrated in the plastic hinges at the region 
boundaries. 
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Angular velocity at the outer boundaries (simply supported) are: 

d_ 0 W dot 
dt Bp r • tan(<j)) 


Internal Energy Dissipation in plate: 


^dot - 4 ■ M 0 • Wd 0 t' 



External work due to uniform pressure p_upper: 

g 

6 := — (aspect ratio of plate) 

L 

Equating these energies and solving for p_upper: 


MlOt ' 


uui r upjpci 


Pupper - 6 ■ M 0 • 


1 + 




tan 


.(♦); 


Bp r 2 .(3-ptan(<|,)) 


Assumptions are made to solve for <|>, the angle of the plastic hinge line along which the velocity profile 
follows. Namely, the minimum p_upper occurs when the derivative of p_upper with respect to tamf> is 
zero. Solving for tan<|), get: 


tan(<(>) = -p + \3 + p 2 which can then be substituted into the energy solution for p_upper 


6 • M 0 . 

Pupper_ss ~ Pupper ss =2.83lx 10 Pa 

Bpr 2 • (^3+p 2 - p) 

Partially restrained supports (applying moment m to outer bound): 0 < m < M 0 * 

m_supp determines the level of restraint on supports: 

m=0 for simply supported, m=Mo for fully clamped ;= 0.50- M 0 


Gives additional internal energy dissapation around plate boundaries: 


Ddot b = 4 m - W dot • 


Bpr 


l") 


+ — 


Bpr ■ tan(<|>) B J 


Add this to internal energy for the simple support case (from above) for: 


Bdot tot - 4 ■ M 0 • Wdot' 


L 
— + 




B tan(<)>) J 



jn_) 

M 0 J 
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The external work from p remains the same as above. Equating and solving for p: 
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DDG-51 Panel from Frame to Frame and Stiffener to Stiffener 


Plate Characteristics: 

Thickness: H;= 0.3438- in 


L : =80ft 


B:= 27.39- in 


~ 4 w ( x >y) + 2• d 2 d 2 w ( x >y) +w ( x »y) = 

dx rlv Hv Hv D 


dx dy 


t 

dy 4 


For HSS Components: 


Cy := 51 • ksi 


ten : = 70 ■ ksi 

t 3 


D := 


EH 


ay = 3.516x lO^kgm 1 s 2 


E rup 0.2: 


12- (l - v 2 ) ^' 2 '<> + v) 


E := 29.6- 10^ ■ ksi E = 2.041 x 10^ kgm * s 2 


p 7833- — 
3 
m 


v := .3C 




( a y + a ten) Plastic Flow Stress 


Blast Characteristics: 


-3 


p 0 := 1 • 10 • Pa t := 0.2- 10 • s 


i ia 6 i - 1 -2 

p 0 = 1 x 10 kgm s 


Static Plastic Collapse Pressure of Rectangular Plates (Jones p.37) 
Plastic collapse moment of solid cross-section: 


■BH 2 


M 0 := ■ 


Plastic collapse moment of solid x-section of unit width: 


Simple Supports: 


Mk :_ ' 


Jo-H 2 


M 0 = 7.952x 10 kgms‘ 



X 

> 

y 

t 

_ _] 










Lower bound (Johansen Yield Condition): 


Define: 



( 


Mx(x) := M 0 


1 -■ 


V 


hr") 


Governing Equation: 


( 


M y (y) :=M 0 


1 -• 






2 "l 

V 1 


1 hr ) 


Mxy( x >y) := 


-Mp-x-y 
Bpr ' Lpr 


2 2 
5'M X + 2 ■ ——M X y + -4-M + p = C 

dx dxdy dy 


Substituting and solving for p: 


Plower ss 2 • M 0 • 


i L f L V 
i + — + — 

B \B) 


5 -1 -2 

Plower ss= l-796x 10 kgm s 


Upper bound (using energy method and kinematically admissable transverse velocity profile): 


Velocity Profile: 


Region I (triangles at ends) Wdot = W dot • 


Region II (center section) 


wdot = W d0 f 


(Bpr • tan(<|>) - Xp r ) 
Bp r • tan 4) 

(Bpr - y) 


B, 


P r 


Each of these regions remains rigid with plastic flow concentrated in the plastic hinges at the region 
boundaries. 
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Angular velocity at the outer boundaries (simply supported) are: 

d_ e w dot 
dt Bp r ■ tan((j)) 


Internal Energy Dissipation in plate: 


Ddot - 4 • M 0 ■ W dot • 



External work due to uniform pressure p_upper: 

B 

p := — (aspect ratio of plate) 

Equating these energies and solving for p_upper: 


Pupper = 6 • M 0 • 


1 + 


p ^ 


tan 




Bp r 2 (3-P tanfo)) 


Assumptions are made to solve for <j>, the angle of the plastic hinge line along which the velocity profile 
follows. Namely, the minimum p_upper occurs when the derivative of p_upper with respect to tan<|> is 
zero. Solving for tan<|>, get: 



which can then be substituted into the energy solution for p_upper 


Pupper ss 



6-M 0 

G/3 + p 2 -p) 


Pupper_ss ~~ l-825x 10 Pa 


Partially restrained supports (applying moment m to outer bound): 0 < m < M 0 " 


m_supp determines the level of restraint on supports: 
m=0 for simply supported, m=Mo for fully clamped 


^upp 0*5 • M 0 


Gives additional internal energy dissapation around plate boundaries: 


f 


Ddot b = 4 • m ■ W d 0 t • 


Bpr 




+ — 


v Bpr • tan (<t>) B J 
Add this to internal energy for the simple support case (from above) for: 


^dot tot - 4 ■ M 0 • Wdot' 


i + _L_V 

B tan (((>) J 


m 0 ; 


in 


0 





The external work from p remains the same as above. Equating and solving for p: 


Ppc 6 ■ Mo 


■^upp ^ 


Bp r 2 - W3+p 2 -p 


p pc = 2.737x 10 5 kgm' 1 s' 2 


For the fully clamped case, the lower bound is: 


Plower fc4 • M 0 ■ 


5 -1 -2 

Plower fc = 2.843x 10 kgm s 


For the fully clamped case, the upper bound is: 


Pupper fc : ~ 


12- M 0 


Bp r ^ • V\ 3 + p 2 - P 


5 -1-2 

Pupper fc = 3.649x 10 kgm s 



DDfXl Panel from Frame to Frame and Deck Gitdet to Deck Girder 


Plate Characteristics: 

Thickness: j^ : = 0.4375- in 


L:=7-54-ft 


B := 9.99 


j4 .2 ,2 A _ n 

—4 w(x,y) + 2 • — 5 —2 w ( x ’ y) + w ( x >y) = — 

dx dx dy dy D 


For HSS Components: 


0y :=51 • ksi 


cr t en : = 70 ' ksi 
.3 


cry = 3.516x lO^kgm ' s 2 


£ rup 0.2; 


D :=- 


EH 


12 -(l - v 2 ) 




2 ■ (l + v) 


E:= 29.6- 10 - ksi 
kg 


p := 7833 


(dy + dten) 


v := 


Blast Characteristics: 

Po := 1 - 10 6 • Pa t := 0.2-10 _ 3 • s 


. ,n 6 i 'I -2 

Po = 1 x 10 kgm s 


Static Plastic Collapse Pressure of Rectangular Plates (Jones p.37) 
Plastic collapse moment of solid cross-section: 


M 0 :=- 


c 0 • B ■ H 2 


Plastic collapse moment of solid x-section of unit width: 


Simple Supports: 




Jo-H 2 



X 


7 

y 

t 




11 - 1-2 

: 2.041 x 10 kgm s 

.3C 

Plastic Flow Stress 


M 0 = 1.288x lO^kgms 







Lower bound (Johansen Yield Condition): 


Define: 





* 




f 


Mx(x) := M 0 


1 - 


V 



Governing Equation: 



M y (y) := M 0 


2 ^ 


1 -• 




M xy (x,y) 


2 2 
-^M x + 2 • M xy + -^,M y + p = C 
0x SxSy dy 


-Mp-x-y 
%r ‘ Lpr 


Substituting and solving for p: 


Plower ss ^ • M 0 ■ 


i L f L \ 
B U; 


Lpr 


4 -1 -2 

Plower ss = 4-534x 10 kgm s 


Upper bound (using energy method and kinematically admissable transverse velocity profile): 


Velocity Profile: 


Region I (triangles at ends) w dot = W dot 
Region II (center section) w dot = w dot 


(Bpr • tan 4) - Xp r ) 
Bp r • tan(([)) 

(Bpr - y) 


B, 


’Pr 


Each of these regions remains rigid with plastic flow concentrated in the plastic hinges at the region 
boundaries. 
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Angular velocity at the outer boundaries (simply supported) are: 

d _ 0 _ w dot 


dt Bp r • tan((|)) 

Internal Energy Dissipation in plate: 


^dot ~ 4 ' M 0 * W^ot g + co ^( ( t ) )j 


External work due to uniform pressure p_upper: Ed ot = 2 • B pr 2 • W dot • p upp er • {- - 

\ B 3 ) 


B 

R := — (aspect ratio of plate) 

T 


1 + — 

Equating these energies and solving for p_upper: p U pper = 6 • M 0 ---— 

Bp r 2 ■ (3 - p 


_P_JI 
.(♦); 


Bp r 2 ■ (3 - p • tan (((>)) 

Assumptions are made to solve for <|>, the angle of the plastic hinge line along which the velocity profile 
follows. Namely, the minimum pjjpper occurs when the derivative of p_upper with respect to tan<t> is 
zero. Solving for tan<|), get: 


tan(<|>) = -p + -J 3 + p 2 
6-M 0 


Pupper_ss 


Bp^-Qs+p 2 -?) 


which can then be substituted into the energy solution for pjjpper 


4 

Pupper_ss = 4.552x 10 Pa 


Partially restrained supports (applying moment m to outer bound): 0 < m < M 0 " 


m_supp determines the level of restraint on supports: 
m=0 for simply supported, m=Mo for fully clamped 


n^upp ^-5 * M 0 


Gives additional internal energy dissapation around plate boundaries: 

Bpr Ll 


Ddot_b — 4 • m ■ W dot * 


Bpr • tan(<|)) B J 


\ r / 

Add this to internal energy for the simple support case (from above) for: 

D do t tot = 4- M 0 - W dot • f-t + -\-r] • f 1 + — 

V B tan4)J 


( 

1 

V 


m 0 ; 
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The external work from p remains the same as above. Equating and solving for p: 


Ppc 


:= 6 • M n 


^ x + m supp ^ 
v. M o ) 


Bpr 2 ' Q3+(3 2 - p) 


4 -1 -2 

p pc = 6.828x 10 kgm s 


For the fully clamped case, the lower bound is: 

M 


Plower fc4 • M 0 • 


4 -1 -2 

Plower fc = 6.123x 10 kgm s 


Bpr 


For the fully clamped case, the upper bound is: 

12-Mo 


Pupper fc 


Bp r 2 • (^3 + p 2 - p) 


4 -1 -2 

Pupper fc ~ 9.104x 10 kgm s 
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DD(X1 Panel from Frame to Ftame and Stiffener to Stiffener 


Plate Characteristics: 

Thickness: H:= 0.4375- in L:=7.54-fi 

A/\Af AAA 


B:= 23.99- in 


.4 ,2 ,2 ,4 

—2 w(x,y) + 2 • — 2 — 2 w ( x >y) + — 4 w(x,y) = — 
dx dx dy dy D 


For HSS Components: 


Gy := 51 • ksi 


G te n : = 70 ■ ksi 
T 3 


8 , -1 -2 


G y = 3.516x 10 kgm s 


£ rup •“ 


E:= 29.6- 10 3 - ksi E = 2.041 x 10 H kgm' 1 s' 2 


p := 7833 • — 
3 
m 


v := .3C 


D := ■ 


EH 


(,_ v 2 ) a= 2(l + ») 


12 


Go := 


(w y + a ten ) Plastic Flow Stress 


Blast Characteristics: 


„-3 


. 6 , -1 -2 


p 0 := 1 • 10 • Pa x := 0.2 ■ 10 • s p 0 = 1 x 10 kgm s 


Static Plastic Collapse Pressure of Rectangular Plates (Jones p.37) 
Plastic collapse moment of solid cross-section: 


M n 


o-B-H 2 


Plastic collapse moment of solid x-section of unit width: 


Simple Supports: 




,-H 2 



x ) 

> 

y 




M 0 = 1.288x lO^kgms 2 


i i. 1 l i TI.I.x i. J 


& 


f 
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Angular velocity at the outer boundaries (simply supported) are: 

d_ Q _ W dot 

dt Bp r • tan(<|)) 


Internal Energy Dissipation in plate: 


D dot = 4 • M 0 * W dot • + cot(<|>) j 


External work due to uniform pressure p_upper: E dot = 2 ■ B pr 2 • W dot - Pupper • ^ 


r : = — (aspect ratio of plate) 
L 


Equating these energies and solving for pjjpper: p upper = 6 ■ M 0 - ■ 


\ tan ((|>)) 

Bpr • (3 - P • tan (cj>)) 


Assumptions are made to solve for 4), the angle of the plastic hinge line along which the velocity profile 
follows. Namely, the minimum pjjpper occurs when the derivative of pjjpper with respect to tan<|> is 
zero. Solving for tanf get: 


a(*) = -P+V3+p 


which can then be substituted into the energy solution for pjjpper 


Pupper_ss ' 


Pupper_ss — 3.764x 10 Pa 


Bpr 2 W3+p 2 -p, 


Partially restrained supports (applying moment m to outer bound): . 0 < m < M 0 ' 


m_supp determines the level of restraint on supports: 
m=0 for simply supported, m=Mo for fully clamped 


®supp 0-5 ■ M 0 


Gives additional internal energy dissapation around plate boundaries: 
r». . _ ^ „ w ' f ^ 


Ddot b = 4 • m • w dot • 


Bpr • tan(<|>) B J 


Add this to internal energy for the simple support case (from above) for: 
Ddotjot = 4 ■ M 0 • W dot • 1 • f 1 + 

V B tanwj V M o ) 
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The external work from p remains the same as above. Equating and solving for p: 


\ + msu PP ^ 


Ppc •- 6 ■ M 0 


V iyi <> ) 


Mo 


Bp r 2 ■ C/3+p 2 - p) 


Pp C = 5.646x lO^kgm ' s 2 


For the fully clamped case, the lower bound is: 

U 2 ) 


Plower fc4 • M 0 • 


5 -1 -2 

Plower fc = 5.939x 10 kgm s 


Bpr 


For the fully clamped case, the upper bound is: 

12 M 0 


Pupperfc ' 


Bp r 2 • a 3 +p~ ~ p 


5 -1 -2 

Pupper fc = 7.528x 10 kgm s 
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APPENDIX G: CALCULATIONS BASED ON DDS100-7 

Calculation for Panel from Frame to Frame and Deck Girder to Deck Girder 

‘ Uniform load over a x b plate - values at center : 

2 _ 4 

Reacts ov = y-Pb 


P-P-b -a-Pb 

a max“ 9 ^max“ ^ AVW ^ L max" 

h Eh 

Maximum stress Maximum displacement Maximum reaction at long 

side (what stiffener must 
support) 

Using flow stress for amax and solving for P: 

K h 

— = 0.876 <—These get |3 —> ~ = 0.755 

S toKlA 


p g := .2874 


(from table) 

P x := .2874 


• 

Must redefine sizes in English units without unit calculations 


hg := .3438 inches 

h x := .4375 inches 


bg := 109.56 inches 

b x := 119.88 inches 


lg 96 inches 

l x := 90.48 inches 

• 

c flow := 6.0510 4 

psi 


a flow b g 

p " nif6 “ »g.p. 

°flow'^x 

P unifx‘“ i n 

Mx 

• 

P unifg - 227l « 

P unifx = 3361 ^ 


These pressures are 

both in psi 


P dds_g = P unifg'P s * 

P dds_x = P unifx P s * 

• 

P dds_g = 1566x lo6pa 

P dds_x = 2 317x lo6pa 
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APPENDIX H: CALCULATIONS BASED ON DDS100-9 


Ambient Pressure 
P a := 14.7-psi 


Peak Incident Overpressure 
Pj := 130.316psi 


Duration of Positive Phase 
- 3 

t := 0 . 2-10 -sec 


Max Reflected Overpressure Dynamic Pressure Stagnation Pressure Peak Avg Pressure 


V= 2P r 


7 p a + 4 ' p l ^ 
7 P a +P l ) 


5 P 1 


H 7P a +P U 


p s . p i + q 0 P c • P 1 ^. 4 -q 0 


P r = 697 . 538 psi 


q 0 = 182 . 044 psi 


P s = 312 . 36 psi P = 57 . 498 psi 


Get plastic section modulus at midspan, Zpm: 


A totg : T g' B g + A g 


A pltg T g B g 


A totx : 4 x B x + A > A pltx :_ 4 x' B > 


= 5 . 613 in 2 


Ap^ = 9.4 Win 2 


totx 2 

- = 6 . 303 in 

2 


A p,tx= 10 . 496 m 2 


Both cases have Atot/2 < Aplt so they correspond to regime 1 in Table 3: 


-namg- 2 . B „ 


7 = 1 r a 2 

mmg- 2 B g' a namg 


7 _ r> j 

lmx’ 0 0 x u namx 


j 2 

Z 2 mg = T' B g'( T g _ d namg) + A webg' 


D g" TF g ' tf ^ 

2 +T g d namgj + A flag'(^ D g 2 + T g d namg j 


1 n (rr j a2 [Dx-^x _ 1 (_ TF x ^ 

Z 2mx- = 2 \ x _ d namv + A webx| 2 + T x “ d namx) + A flax D x + T x “ d namxj 


z PMg Z lmg + Z 2 mg 


^PMx ■ ^ 1 mx + ^2m^ 


z PMg = 8 - 7 45 in 3 


Z PMx = 12 . 874 in 
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Natural period of stiffeners 

M tg - = ( B g' T g + A g)‘ L g’ p ms M tx := ( B x' T x + A x) ,L x'P 


ms 


M t = 304.991b 


M tx = 323.0171b 


K ebarg 


307-Ejns-Ig 


TR 




K ebarx ■ 


w 


for \x > 1 


k ebarg = 2 - 578 * «> 7 

sec 


^ebarx 6.71 x 10 


7 lb 


sec 


^lm 0.6f for jit > 1 


T ng 2 ' n 


f K LM" M tg ^ 

kebarg J 


0.5 


T nx := 2 : n ' 


^ ^ebarx ) 


0.5 


T n g = 0.018sec 


T nx = 0.01 lsec 


For single pulse loading case, required resistance: 


P _ T ng/_, ,\°-5 2-^1-1 

F ig : ~~'( Fl ) + - 

2 Hj- 


1 + 

( T VI 
0.7--^ 

r lx-- 

n-T 


1 x )\ 



( 2 '“1 - >) 


0.5 


2(ij - 1 




( T Yl 

*nx 1 

2 Hj- 

1 + 

0.7— 


_ 

l t Jj 


F2g:=^.( 2 M2-r + 

& n .-[ V 1 


0.5 


2-n 2 - 1 


2-h 2 ' 


1 + 

( T Yl 

0.7--^ 


1 t J. 


F 2 x :=-( 2 . 2 - 1)° 5 + 

7T*T x 7 

P s' B x^ 


“rig 5 " 




P s B g L g 

F lg 

P s B g- L g 

F 2g 


^g 


= 2.49 x 10 newton 


Rj. lg = 6.425x 10 lbf 


R,-2 g = 1.455x 10 lbf 


Rj.|g = 2.858x 10 newton 


2-h 2 - 1 



f T nxYl 

2-Hj 

rK 


R r lx := 
R r2x := 


F lx 

Ps-Bx 1 ^ 

F 2x 


Rj.j x = 8.319x 10 lbf 

Rj. 2x = 1.885x 10 4 lbf 
„5 


5 Rj.| x =3.7x 10 newton 

Rp X = 3.702x 10 newton 

^r2g = 6*473 x 10 4 newton Rj. 2 x = 8.383 x 10 4 newton 



APPENDIX I: COMPARISON OF STIFFENED PANEL 

CALCULATIONS 








DDG - 51 Class 

DD(X) 



in 

0.3438 

0.4375 


L (Frame Spacing) 

ft 

8 

7.54 



ft 

9.13 

9.99 


Stiffeners per panel 


3 

4 



ft 

3.043 

2.498 


Static Failure Pressure (Moment 




W 

calculation) 

Pa 

1.23E+04 

1.67E+04 

0) 

E 

Dynamic Failure Pressure (Moment + 




E 

strain rate effects) 

Pa 

3.65E+04 

4.93E+04 

£ 

Simply Supported Lower Bound 




0) 

■o 

(Johansen Yield Condition) 

Pa 

2.83E+04 

4.53E+04 

O 

Simply Supported Upper Bound (Energy 




■o 

Method with Velocity Profile) 

Pa 

2.83E+04 

4.55E+04 

c 

re 

Fully Clamped Lower Bound (Statically 




V) 

JC 

Admissible Stress Field) 

Pa 

3.78E+04 

6.12E+04 

u 

0) 

Fully Clamped Upper Bound (Energy 




Q 

c 

Dissipation Method) 

Pa 

5.66E+04 

9.10E+04 

0) 

re 

Partially Clamped (50%) (Energy 




2 

Method) 

Pa 

4.25E+04 

6.83E+04 

m 


Pa 

1.57E+06 

2.32E+06 

IS 

DDS 100-9 (air blast / nuclear) 




CL 

For Max Deflection w/o rupture 

Pa 

2.55E+07 

2.94E+07 


For Ship Specifications 

Pa 

5.75E+06 

6.66E+06 


Static Failure Pressure (Moment 




<S 

p 

calculation) 

Pa 

4.93E+04 

1.04E+05 

c 

E 

Dynamic Failure Pressure (Moment + 





strain rate effects) 

Pa 

1.46E+05 

3.08E+05 

E 

<D 

Simply Supported Lower Bound 




*D 

u 

(Johansen Yield Condition) 

Pa 

1.80E+05 

3.71 E+05 

CD 

~a 

Simply Supported Upper Bound (Energy 




c 

re 

Method with Velocity Profile) 

Pa 

1.83E+05 

3.76E+05 

E 

Fully Clamped Lower Bound (Statically 




o> 

c 

Admissible Stress Field) 

Pa 

2.84E+05 

5.94E+05 

I 

Fully Clamped Upper Bound (Energy 




w 

Dissipation Method) 

Pa 

3.65E+05 

7.53E+05 

c 

<u 

Partially Clamped (50%) (Energy 




1 

Method) 

Pa 

2.74E+05 

5.65E+05 

re 

m 


Pa 

2.48E+06 

4.53E+06 


DDS 100-9 (air blast / nuclear) 




re 

0L 

Ml 1 11 UWIUdMUllMIIU 

Pa 

2.55E+07 

2.94E+07 


For Ship Specifications 

Pa 

5.75E+06 

6.66E+06 
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